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1. Introduction and General Status 

During the fourth semi-annual period (12/30/89 - 6/30/90) of this project our work 
has concentrated on the solution of a two-dimensional axisymmetric moving boundary 
problem for the sensitivity of isothermal and non-isothermal liquid columns and the 
sensitivity of thermo-capillary flows to buoyancy-driven convection caused by residual 
accelerations. The tasks proposed for this period are outlined in our second year 
continuation proposal. A summary of progress to date is given below. 

2. Progress Summary 

The following tasks were described in the original and continuation proposals. 

Task 2 : Thermo-capillary Convection - Two Dimensional (fixed planar surface) 

Task 3: Thermo-capillary Convection - Three Dimensional and Axisymmetric 

Task 4: Liquid Bridge/ Floating-Zone Sensitivity 

Task 6: Transport in Closed Containers 

Task 7 : Interaction: Design and Development Stages 

Task 8: Interaction: Testing Right Hardware 

Task 9: Reporting 


Work proceeded roughly according to schedule. A review of work related to 
experiment sensitivity to residual acceleration was completed and has been accepted for 
publication. A copy of the article is attached in Appendix 1. Much of the 4th semi-annual 
period was spent completing computer codes for the solution of the liquid bridge free 
surface sensitivity problem and accumulating results for Task 2. (This was done in 
anticipation of the hiatus in funding which prompted the extension of the period of 
performance for the second year’s work.) In the third semi-annual report we stated that 
Task 2 should complete by June-July of 1990 (i.e. earlier than originally anticipated). Since 
that report was issued NASA MSFC has extended the P.O.P. of the second year’s work to 
June 1991 as continuation funds are not yet available. We now anticipate completion of 
Task 2 during this period. Task 2 involves an analysis of the residual acceleration 
sensitivity of thermo-capillary driven flow in a 2-D rectangular region. The free surface is 
held planar. Progress to date will be briefly summarized in Section 4. 

Task 3, which involves the examination of axisymmetric thermo-capillary 
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convection has been started early in the process of developing a method for the examination 
of non-isothermal liquid bridge. A code for the axisymmetric isothermal and non- 
isothermal (i.e. thermo-capillary driven convection is included along with a deformable free 
surface) liquid bridge sensitivity analysis has been developed and a comparison with 1-D 
results has been made. These results will be reported in Section 3. 

Two papers summarizing work carried out under the VIT project have been 
accepted for publication during the 4th semi-annual period: 

1. The Sensitivity of a Liquid Bridge to Axial Vibration, with Y. Q. Zhang, 
accepted for publication, Physics of Fluids, 1990. 

2. Microgravity Experiment Sensitivity to Residual Accelerations: A Review, 
accepted for publication. Microgravity Science and Technology, 1990. 

One paper is near completion: 

3 . A Finite Difference Method for a Model Float Zone Free Surface Problem, 
with Y. Q. Zhang, to be submitted to the Journal of Numerical Methods in 
Fluids, 1990. 

A fourth article concerning time dependent response of isothermal and non-isothermal liquid 
bridges is planned. 

3. Liquid Bridge/Floating-Zone Sensitivity 

This section refers to Task 4 and concerns the sensitivity of an isothermal liquid 
bridge to axial acceleration and has been extended to non-isothermal situations covered 
under Task 3. The first annual report contained a detailed description of a one-dimensional 
model and the sensitivity results obtained using that model, this work has subsequently 
been refined and will shortly be published in the journal Physics of Fluids. A copy of the 
accepted manuscript is attached in Appendix 2. 

The 1-D model had been used to examine the sensitivity of an isothermal liquid 
bridge to axial vibration. The results indicate that the zone is most sensitive to accelerations 
with frequencies close to or equal to the lowest natural frequency of the zone. For the 
purposes of this project it is most useful to assess the sensitivity of the zone in terms of 
predicted Space Station and/or Spacelab environments. For the cases examined, we have 
seen that frequencies around the 0.1 Hz range appear to be the most sensitive. The low 
frequency (< 10‘ 2 Hz) acceleration environment predicted for the space station 1 should not 
exceed levels of 10'^ g. Higher frequencies can be associated with acceleration magnitudes 
of up to 10' 2 g. In terms of these predicted levels, or those measured on past missions 1 , 
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the practical sensitivity range is restricted to disturbances with frequencies ranging from 
10- 1 - 10 Hz. 

During the second year, work on a full axisymmetric calculation started. A model 
was developed for both isothermal and non- isothermal bridges. The adaptation of two 
numerical algorithms for the solutions of these problems has been successfully completed 
and preliminary calculations have been made. 

The early stages of development and verification of the isothermal code (and the 
ability of the method to handle buoyancy-driven flows in fixed geometries) was outlined in 
the third semi-annual report and an updated account is given below. 

For the solution of the axisymmetric isothermal and non-isothermal problems our 
objectives were: 

1 . To develop a full numerical model of the axisymmetric free surface 
response of an isothermal liquid bridge (suspended between rigid disks of 
equal radius) subject to axial vibration and to compare it to the sensitivities 
predicted by a previously developed 1-D model 

2. To develop a full numerical model of the axisymmetric free surface 
response of a non-isothermal liquid bridge subject to axial vibration. The 
presence of a temperature gradient gives rise to a thermo-capillary flow. We 
wish to investigate the interaction between the free surface motion induced 
by vibration, and the surface-driven flow. Interaction of these flows with 
internal buoyancy-driven flow driven by the vibration is expected to be 
negligible but is included in the basic model. 

Any algorithm that is capable of solving the equations governing the second 
objective must be capable of solving the first objective. Added complications result from 
the fact that the coupled heat transfer and Navier-Stokes-Boussinesq equations must also be 
solved to satisfy the second objective. In addition, it is convenient to use steady-state 
solutions for thermo-capillary convection as initial conditions. As a result we have 
developed two algorithms which are based on the same numerical methods; one is used to 
compute steady-state velocity and temperature fields and free surface shape for the 
following problem, the other is used to calculate the time dependent response. 
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3.1 Steady State interface shapes and thermo-capillary convection 

Consider a cylindrical liquid zone (see Fig. 1) held between two parallel coaxial 
circular rigid disks (radius = Ro) separated by a distance L. The liquid is a non-isothermal 
incompressible Newtonian fluid. The bridge is held between the disks by surface tension. 
The free surface of the bridge is a gas-liquid interface and is described by r=R(z,t). The 
two disks are maintained at a constant temperature T 0 . Surface heating is provided through 
a parabolic ambient temperature Too(z) and the heat transfer coefficient at the free surface is 
denoted by h. In addition, we make the assumptions that the gravitational acceleration is 
parallel to the cylinder axis and that the velocity, temperature and free surface deformation 
are axisymmetric. Furthermore, we let the surface tension at the free surface vary linearly 
with temperature and assume that the Boussinesq approximation holds. 

The governing equations are made dimensionless by scaling length, time and 
velocity with Rq, Rq/U* and U*> respectively . Here U* is a characteristic velocity given by 

„• l?|AT 
H 

where AT=T max -T m i n represents the maximum temperature difference at the surface, lyl is 
the absolute value of the derivative of the surface tension with respect to temperature, and p 
is the dynamic viscosity. We shall refer to a "half zone” model when the ambient 
temperature has extrema at the end disks and a "full zone" when the temperature maximum 
occurs between the disks. For a full zone we shall take T max to be T<w(0), and T m i n to be 
T 0 , where A = L/Ro is the aspect ratio. 

The non-dimensional pressure is 

_ P* + Pogz D 

p = — — Ro 

poU 2 

where p* is the dimensional pressure, g is the gravitational acceleration, z is the 
dimensionless axial coordinate, and po is the density corresponding to the reference 
temperature. The temperature is rendered dimensionless using T max - T m i n - With these 
scales the dimensionless equations in a cylindrical coordinate system can be written as 
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where the Reynolds number. Re, Marangoni number, Ma and Grashof number Gr, are 
respectively 


Re-B*Sl . M.-UISJ5.. Gr-iS"*. 

V JJ.K V 2 

Here, v is the kinematic viscosity, K is the thermal diffusivity, P is the volume thermal 
expansion coefficient, and g is the gravitational acceleration. At the disks, the boundary 
conditions are 

u = w = T = 0, at z = ± (5) 


The symmetry conditions at the centerline r = 0 are 

dw dT 

U “ ~dF~ dF~‘ 


( 6 ) 


The free surface is located at r = R(z) where the boundary conditions are: 
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and are, respectively, the capillary number, Biot number, and dimensionless gravitational 
acceleration, y 0 is the mean surface tension, and k is thermal conductivity. The force 
balance at the free surface in the normal and tangential directions are given by eqs. (7) and 
(8), respectively. Equation (9) is the kinematic boundary condition at the liquid-gas 
interface. The thermal boundary condition at the interface is given by equation (10) in 
which the equivalent heat transfer coefficient h contains the effect of the radiant and 
convective heat transfer. The constant X in (7) represents a dimensionless reference 
pressure difference 2 ’^ across the interface. In model float zone systems with fixed rigid 
endwalls, such as the one discussed here, X is determined by the following constant 
volume constraint 
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Finally, the condition that the contact lines between the liquid end disks are fixed is 


R = 1 at z = ± A/2. (12) 

Numerical Method 

In the case of a two dimensional axisymmetric flow, the governing equations can be 
simplified by introducing the stream function, and vorticity, to, as new dependent 
variables: 
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Substitution of (13) into (14) yields 


rco = 


3V i dy 3V 
a7rar + a?- 


(16) 


The original set of three equations governing mass and momentum has been reduced to 
two equations governing the steam function and vorticity. 

The steady free boundary problem for cylindrical liquid zone is solved iteratively, 
since the location of the free surface is a priori unknown. To obtain a solution we adopt a 
Picard iterative procedure 4 as follows: 

1. Guess the free surface shape for the initial iterate; 

2. Obtain the approximate temperature and velocity fields by transforming the governing 
equations and boundary conditions to a circular cylindrical domain via a non-orthogonal 
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transformation and solve them using a pseudo-unsteady semi-implicit method; 

3. Obtain the pressure at the free surface by integrating the transformed momentum 
equation; 

4. Use the normal force balance condition at the free surface to decide how to update the 
free surface location; 

5. Return to step 2. Repeat until convergence is obtained. 

N on-orthogonal Transformation 

The region occupied by the liquid zone is transformed into a fixed circular 
cylindrical computational region using a non-orthogonal coordinate transformation, i.e. 


T] = Z , 



(17) 


It then follows that 


3 i 3 3 3 C 3R d 

dr R ^ ’ dz 3rj R 3q ^ 


The advantage of this transformation is that the free boundary coincides exactly with a 
coordinate line in a computational grid and regeneration of mesh during the outer iteration 
is avoided. For the sake of simplicity, a pseudo-unsteady method in association with a 
semi-implicit discretization with respect to time is used. 

The resulting system of governing equations is solved as the time (pseudo) 
derivatives of T, y and co — » 0. The discretization in time employs an explicit Adams- 
Bashforth scheme for the nonlinear convective terms and the implicit scheme for the 
viscous terms. Other terms in the equation are treated explicitly. For spatial discretization, 
central differences are used. The A.D.I. form for the vorticity equation is then 
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Here the superscript n+1/2 denotes the intermediate step associted with the ADI 
method 7 . The velocities (at the nth step) are taken from the values at the previous step 
values. Thus, except for the thermal condition at free boundary, these are Dirichlet 
boundary conditions. The heat transfer equation is also discretized in an ADI form. The 
resulting system of discretized equations is solved using a factorization method 5 . 

We define a steady state to occur reached when the residuals (dcii/dx, dT/dx, dy/dx) 
of the voracity, energy and stream function equations are less than 10‘ 7 . That is, 


max 


F " +1 . F " 
r U ij 


Ax 


< 10‘ 7 . 


( 21 ) 


Where F represents the vorticity, temperature and stream function at the n+lth and nth 
iterative step. The numerical solution is second order accurate in space. 

Having computed the vorticity and stream functions for a given surface shape, it 
remains to iterate on the condition for the balance of force normal to the surface in order to 
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obtain the final steady surface shape. In addition, the shape must satisfy the volume 
constraint (11) and boundary conditions. 

Two iterative schemes for determining the interface are discussed in this section. In 
both schemes, the new velocities and temperatures are taken from the current calculated 
values and the pressure at interface can be obtained by taking the r| component of 
momentum equation and integrating with respect to r\. This yields 
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The first scheme (I) involves successive approximation by the direct solution of 
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This scheme is based on the following principle. A shape is assigned to the free surface with the 
calculated pressure, velocity and temperature. An initial guess for the pressure constant X is made. 
The new interface shape is then determined directly from (22). The integral (1 1) is then evaluated 
to check whether the volume constraint has been satisfied. If it is not satisfied an inner iteration is 
made using a Newton-Raphson procedure to calculate the following improved estimate of X.: 

X k+1 = X k - (— V*AV , (23) 


where 
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M/2 

AV = I 7tR 2 dr| - V 0 . (24) 

J- A/2 

The above procedure for determining X is quite effective and is repeated until the volume 
constraint is satisfied. R(rj) is then updated. New velocity, pressure and temperature fields 
are calculated using the updated value of R(r|). The outer iteration is repeated until 

max |R n+1 - R"! < e , (25) 

where we take e = 10* 4 . 

The second scheme (II), used by Ryskin and Leal 7 , uses the residual of the normal 
component of force balance to drive the shape to its steady position. This is equivalent to 
equating the residual with an artificial capillary force. This effective force causes a local 
displacement of the interface in the direction of the force. The magnitude of the local 
displacement is proportional to this force. The interface shape at each iteration is thus 
modified so as to reduce the residual until the condition (22) is met. It follows that at each 
iteration the improved interface shape is given by 

Rj m+1 = R m + ocExj, (26) 

where Exj is the residual of the force balance equation at the jth surface location and the 
constant coefficient a is determined by numerical experiment. In order to ensure 
convergence, a should be small. We found that the values of a which led to rapid 
convergence depended on the product of Re* 1 and Co' 1 . If a is chosen to be too small the 
amount of CPU time used increases substantially. 

The change in volume between the mth and (m+l)th iteration can be found from the 
volume constraint (11) and equation (26) and neglecting higher order terms, i.e. 

M/2 

I Exj Rj 11 dr| = 0. (27) 

J A/2 

The pressure constant X is contained in Exj and is obtained by satisfying (24). Even then 
the liquid bridge may still change volume slightly at each iteration owing to numerical error 
and higher-order effects. These small changes can accumulate and eventually result in a 
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gross error. To prevent this, the formula (26) is modified to 

Rj m+1 = Rf j^) 1/2 + aExj . (28) 

We have used the method described earlier to examine the influence of various 
parameters on momentum and heat transport and meniscus shape. In addition, a 
comparison between the results obtained with the combined scheme and those of Hyer et 
al. 7 has been made. 

A parametric study of the effects of varying the temperature difference AT (keeping 
Ma/Gr fixed), Marangoni number, Ma, Nusselt number, Nu, aspect ratio. A, and Grashof 
number, Gr, has been made and will be described in a forthcoming paper. 

Unsteady interface shapes (g-jitter response) for isothermal and non- 
isothermal bridges. 

To obtain the time dependent internal and free surface response of isothermal and 
non-isothermal liquid bridges to axial acceleration we employed the following procedure. 

To obtain a solution we follow the approach of Kang and Leal 8 and adopt a Picard iterative 
procedure 4 as follows: 

1. Use the free surface shape, stream function, vorticity, velocities, and temperature 
computed from the steady calculation as the initial condition. (For the isothermal case the 
liquid bridge is initially static with a free surface shape given by the prevailing steady axial 
acceleration ; 

2. To obtain the temperature and velocity fields at the next time step we transform the 
governing equations and boundary conditions to a circular cylindrical domain via a non- 
orthogonal transformation and solve them using a semi-implicit Adams-Bashforth/Crank- 
Nicolson time discretization. Centered differences are used for the spatial approximation ; 

3. Obtain the pressure at the free surface by integrating the transformed momentum 
equation; 

4. Use the normal force balance condition at free surface to determine the updated free 
surface location at the new time step; 

5. Return to step 2. Repeat until convergence is obtained. 


Our early results (discussed in the 3rd semi-annual report) showed good agreement 
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between the predictions of the 1-D isothermal model and the axisymmetric calculations (see 
Table 1) Further calculations have been made and for the cases examined to date these 
findings are confirmed. For the cases examined to date the effect of thermo-capillary 
convection on the response of the free surface shape has been found to be negligible (See 
Fig. 2.) Figure 3. depicts sensitiviies of non-isothermal bridges for A = 2 and 2.6. 

Table 1 : Comparison of Bridge Radii (R) (to 3 decimal places) at two locations predicted by 
1-D and axisymmetric (2D) models, g=l. 42x10-3 ms -2 , the frequency of axial vibration is 
3 Hz. 


Time (s) 

0 

0.22 

0.39 

0.48 

Ri-d 

1.001 

1.149 

1.017 

0.891 

ESBBM 

1.001 

1.146 

0.998 

0.919 

Ri-d 

1.001 

1.159 

0.994 

0.909 

Raxisym 

1.001 

1.150 

0.986 

0.929 
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FREQUENCY (hz) 


Figure 3. The acceleration sensitivity predicted by the axisymmetric non-isothermal 
model for a liquid bridge with A = 2 (solid line) and A=2.6 (dashed line), with Ma - 

39.75. The liquid has the properties of molten silicon. 
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4. Thermocapillary Convection -2D 

The work for Task 2 involves the solution of the incompressible Navier-Stokes - 
Boussinesq equations in a rectangular region with one "free surface" boundary condition 
which equates the tangential viscous force at the surface to the temperature gradient along 
the surface. A pseudo-spectral Chebyshev collocation method^ has been employed. The 
basic approach and the equations to be solved were oudined in the first annual report. 

The objective is to study the effect of g-jitter on thermo-capillary flows for a range of 
dynamic Bond numbers Bd = pgL 2 /[cty/9T] (here the denominator represents the change in 
surface tension with respect to temperature). To date we have examined fluids with Prandtl 
numbers in the range 10' 2 to 10, which spans the range associated with most experimental 
materials. Marangoni numbers in the range 1 to 10 3 have been subjected to single 
frequency excitations with amplitudes in the range ltH-lO' 2 g. As expected, the sensitivity 
of the flow to a given residual acceleration depends on the relative magnitudes of the 
Reynolds number and the Grashof number (see discussion in the 1st annual report and 3rd 
semi-annual report). 

Description of the model system (see also 3rd semi-annual report) 

The physical system currently under examination involves a 2-D rectangular region 
which is bounded by three rigid walls and a "free surface". Two of the rigid walls are 
isothermal with temperatures Th and Tc. The other boundaries are adiabatic (no 
perpendicular heat flux) and the "free surface" is constrained to be a straight line. This 
retains the essential character of a free liquid surface while removing the numerical 
problems associated with a deformable one. The forces at the interface may be decomposed 
into components normal and tangent to it. Contributions to the force balance include the 
fluid pressure, the viscous forces, and the capillary forces. For the normal component of 
force balance, the difference between the pressure and viscous forces on either side of 
theinterface are proportional to the interface curvature multiplied by the surface tension y. 
For the tangential force balance at the interface, the difference in the shear stress 
components is balanced by a gradient in surface tension. Such gradients will arise when 
there is a temperature gradient along the interface. The gradient in temperature produces a 
gradient in surface tension that induces shear stresses and, thus, fluid (thermo-capillary) 
motion. In our calculations the residual gravity vector is capable of taking any orientation 
with respect to the boundaries and can be time-dependent. 
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Preliminary results were presented in the 3rd semi-annual report. A more detailed 
account of our results will be written during the extended period of performance and will 
include an analysis of low Prandtl number fluids (Pr <1 ), the effects of 1 Hz disturbances 
and an examination of system size effects. In addition the response to a realistic disturbance 
taken from SL-3 residual acceleration data will be examined and discussed in light of our 
understanding of the response of the system to impulse accelerations. 

5. Thermo-capillary Convection - Three Dimensional and Axisymmetric 

This work was not scheduled to start until the sixth month of the second year 
although a 2-D axisymmetric code has been developed as part of Task 4. This is described 
in some detail in Section 3. 

6. Transport in Closed Containers 

In the original proposal we agreed to make available the results of related work at 
the CMMR not funded through this grant. These results we included in the appendices of 
the 3rd semi-annual report. 

7. Work planned for the Extended Period of Performance 7/1/90 - 6/30/91 

As we stated in Section 1., much of the 4th semi-annual period was spent 
completing computer codes for the solution of the liquid bridge free surface sensitivity 
problem and accumulating results for Task 2. This was done in anticipation of the hiatus in 
funding which prompted the extension of the P.O.P for the second year’s work. 

During the extended P.O.P. we plan to complete the work associated with Task 2 
and complete a Final Report. Obviously, work that was to be completed during the third 
year cannot be completed until funds are made available. A continuation proposal for the 
remaining year’s work will be submitted during the extended period of performance. 
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LOW-GRAVITY EXPERIMENT SENSITIVITY TO 
RESIDUAL ACCELERATION: A REVIEW 

J. Iwan D. Alexander 

Center for Microgravity and Materials Research 
University of Alabama in Huntsville 
Huntsville, Alabama 35899 

Abstract 

In this article, work related to the analysis of experiment sensitivity to residual 
acceleration experienced in low earth orbit spacecraft is reviewed. Most of the work 
discussed concerns heat, mass and momentum transfer in fluid physics and materials 
science experiments. On the basis of our current understanding of experiment sensitivity it 
is concluded that, in general, experimenters should be concerned about the effect of 
residual acceleration and that careful modelling included as part of an experiment program 
will enable optimal use of the limited experiment time available in space. 
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1. Introduction 

It has been recognized for some time that the low-gravity acceleration environment 
associated with a spacecraft in low earth orbit offers an opportunity to study certain 
physical processes which are difficult to investigate under the gravitational acceleration 
experienced at the earth's surface. Our experience with such opportunities has led us to 
realize that the residual acceleration environment on board an orbiting spacecraft is not as 
low or as steady as would be desired for certain classes of experiments . The sources of 
the residual acceleration include [1-5] crew motions, mechanical vibrations (pumps, 
motors, excitations of natural frequencies of spacecraft structures), spacecraft maneuvers 
and basic attitude motion, atmospheric drag and the earth's gravity gradient. The 
accelerations are characterized by temporal variations in both magnitude and orientation. 
Such disturbances will, in particular, affect those experiments susceptible to buoyancy 
induced fluid motion [6]. Indeed, recent studies [7,9] indicate that transient disturbances 
can have undesirable long term effects. The analysis of the sensitivity of any proposed 
flight experiment is necessary for a variety of reasons. These include the need to optimize 
the use of the limited time available for flight experiments, the interpretation of experimental 
results and the determination of tolerable acceleration levels to be used in planning for 
NASA's Space Station [3]. 

Examples of experiments that are conducted in the low-gravity environment include 
critical point studies [10], crystal growth [11-18], and diffusion experiments [19,20]. 
Since these experiments involve density gradients in the fluid phase, they are inherently 
sensitive to the effects of buoyancy-driven fluid motion [21-25]. For certain types of 
systems it is known [7,8,26-37] that the associated physical processes (mass transfer, heat 
transfer, convection, oscillation and distortion of free surfaces) are sensitive to time- 
dependent accelerations. For example, the transfer of mass and heat in fluid systems can 
be significantly affected by oscillatory flows [28-37], The effect of oscillation enhanced 


3 


transport conditions on the local variation in composition during crystal growth is not well 
characterized for most of the systems relevant to crystal growth under microgravity 
conditions. Order of magnitude analyses [7,8,38] suggest that for certain combinations of 
physical properties and growth conditions, oscillations in the residual acceleration may 
adversely affect the mass transport conditions. For instance, the DMOS experiment [39] on 
STS 51-D showed extensive evidence of convective mixing of liquids. The degree of 
mixing is greater than can be attributed to the quasi-steady low-gravity component, but can 
be accounted for when oscillations in effective gravity are included in the description of the 
transport conditions. 

It is possible that in future space laboratories some experiments will require isolation 
from vibration. The design of efficient isolation systems, however, requires determination 
of the tolerable acceleration levels (amplitudes and frequencies) for given experiments. The 
results of sensitivity analyses can be used in the development of vibration isolation systems 
to identify those experiments that will need isolation and by supplying data concerning 
frequency dependence, the effect of transients etc.. The system can then be designed to 
filter out those bandwidths to which the experiment is predicted to respond adversely. 

Low-gravity experiments which involve liquids with free surfaces are most likely to 
be sensitive to vibrations. Experiments of this type include studies of liquid bridges, 
equilibrium shapes of drops and bubbles [40-43] and thermocapillary flow experiments 
[43-46]. Evidence for the oscillation of liquid zone shapes was found on the D-l mission. 
Long liquid bridges were more sensitive to residual gravity, and exhibited random 
oscillations in zone shape [40]. 

It should be recognized that the sensitivity of a given process to the overall low- 
gravity environment will depend on one or more of the following factors: 

a) magnitude and direction of accelerations 


b) system geometry 
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c) boundary conditions (eg. insulated or conducting walls, solute sources and sinks), 
and 

d) physical properties of the participating materials (viscosity, thermal and solute 
diffusivities). 

Any analysis of the sensitivity to oscillatory accelerations and transient disturbances 
characteristic of the anticipated environment aboard Space Station, or any other spacecraft, 
must take these factors into account. In general, a full mathematical characterization may 
not be practical. However, identification of sensitive experiments in combination with a 
judicious choice of simplified numerical models can be used to choose operating parameters 
that will optimize the use of the low-gravity environment. With this in mind, we review 
work related to the effects of steady and time-dependent residual accelerations, and 
associated work concerning transport in oscillating flows. 

In section 2 we describe the components of the low-gravity environment. In section 
3 previous work involving order of magnitude (OM) estimates is discussed. The 
application of OM analyses to acceleration- sensitivity determination is discussed in 3.1. 
The results of simple analytical models are summarized in section 3.2. Work related to 
linear stability in the presence of modulated gravity is briefly discussed in section 3.3. The 
results of direct numerical simulations are reviewed in section 3.4. In section 4., the 
results and utility of residual acceleration analyses are summarized and discussed. 
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2. The microgravity environment 

The low-gravity environment of a spacecraft has been referred to as a "zero gravity" 
environment. This is because any object within the craft is subject to roughly the same 
acceleration due to the earth's gravitational field. Thus, a free object may move along 
approximately the same orbital path as the spacecraft. As a consequence, to an observer in 
the space craft objects appear to behave as if no gravity were present. However, any body 
capable of motion relative to the spacecraft will experience an acceleration relative to the 
mass center of the spacecraft. This relative acceleration arises from several sources which 
include the gravity gradient of the earth, atmospheric drag and attitude motions, as well as 
machinery vibrations and crew activities. 

The relative acceleration that is associated with the earth's gravity gradient is defined 
as follows. The mass center of the spacecraft is subject to a force F(r 0 ) due to the 
gravitational attraction of the earth, where r 0 is the position of the mass center of the craft 
with respect to the mass center of the earth. A particle at a position r within the spacecraft is 
subject to a force F(r). If, as a first approximation, we ignore the interaction between the 
spacecraft and the mass of the particle, and if the distance lr - r 0 1 is small compared to the 
semi-major axis of the spacecraft orbit, the difference between the forces applied at r and r 0 
define the gravity gradient G as follows [2-6] 

F(r) - F(r 0 ) = Gz, z = r - r 0 . (1) 

If the spacecraft maintains a fixed orientation, for example with respect to the sun, 
there is no rotation of the spacecraft frame relative to the geocentric frame. In this type of 
orbit both the magnitude and the orientation of the residual acceleration vector change as a 
function of time and there will be no steady residual acceleration. 

In addition to the gravity gradient acceleration the atmospheric drag force on the 
spacecraft can be significant [2-6], despite the fact that the atmospheric density at typical 
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shuttle altitudes is only on the order of 10' 12 kg nr 3 . Atmospheric drag causes a slow 
inward spiral of the spacecraft. To an observer in the spacecraft the effect of atmospheric 
drag is to produce a relative acceleration that is equal in magnitude but opposite in sign to 
the atmospheric drag force. Table 1 gives atmospheric drag acceleration estimates and 
relative acceleration related to the gravity gradient and centrifugal force in pg per meter 
from the mass center of a spaceraft that is in a circular orbit with a gravity gradient 
stabilized attitude. Xi is parallel to the local vertical and X 2 is perpendicular to the orbital 
plane. For this attitude the acceleration points away from the mass center along Xi and 
toward the mass center along X 2 . 

If the spacecraft does not maintain such a fixed orientation then, in addition to (1), 
centrifugal, Coriolis and Euler accelerations will become apparent in the spacecraft 
reference frame. They arise in conjuction with spacecraft attitude motions and depend on 
the nature of the rotation of the spacecraft frame of reference relative to a fixed geocentric 
frame. If a spacecraft in a quasi-circular orbit continuously rotates relative to the fixed 
geocentric frame such that a given direction in the spacecraft frame is always oriented 
parallel to the craft’s position vector r 0 , a centrifugal acceleration will augment the gravity 
gradient in the direction parallel to the position vector, and cancel the component tangent to 
the flight path. Under these conditions there will be a steady component to the residual 
acceleration vector. In practice this situation arises, for example, in the so-called gravity 
gradient stabilized attitude. 

Euler accelerations arise when the rate of rotation of the spacecraft frame changes as a 
function of time, as it would in a non-circular orbit, or as the semi-major axis of the orbit 
changes as the craft slowly spirals inward due to atmospheric drag. Coriolis accelerations 
will also be apparent but are proportional to the velocity of an object relative to the 
spacecraft. Thus, unless relative velocities are large, the Coriolis effect will be insignificant 
[ 6 ]. 
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In addition to the above, there are a variety of residual acceleration components that 
occur over a broad range of amplitudes and frequencies. The nature and sources of these 
accelerations have been discussed in [1-3]. Figures 1 and 2 are examples of residual 
accelerations measured on orbit. 

3. Analyses of residual acceleration effects 

As discussed in the introduction, many experiments are undertaken in orbiting 
spacecraft in order to reduce or eliminate unwanted effects of buoyancy driven motions in 
fluids. Sources of buoyancy in fluids include discrete interfaces, as they occur, for example 
between immiscible fluids, solid particles and fluids, and bubbles and fluids with different 
densities, other buoyancy sources are internal density gradients that arise as a consequence 
of gradients in temperature and concentration. The nature of the response of a bulk fluid to 
the effective body force associated with a low-gravity environment will depend, in part, on 
the relative orientation of the residual acceleration vector and the density gradient. In 
particular, if the integral of the effective body force pb (where p is the density and b the 
effective acceleration) acting on any closed circuit within the fluid is non-zero, motion will 
immediately ensue. That is, if 

j dc-pb = j dA-curl pb = j dA-[grad pAb + p curl b] *0. (2) 

motion will occur. Note also, that even if there is no density gradient motion will occur if 
j dA-curl pb *0. 


( 3 ) 
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This includes the case of "precessional stirring" [48]. This could occur, for example, when 
a spacecraft orbiting in a radial attitude (i.e. it continuously rotates about an axis 
perpendicular to the orbital plane) underwent an additional rotation about some axis other 
than the basic orbital axis. 

There have been three appraoches to the determination experiment sensitivity to the 
residual acceleration environment. 

a) order of magnitude analyses [38,49-52], 

b) analytical models [5,6,39,53-59], and 

c) direct numerical simulation [9,3,60-73]. 

In addition to work directly related to microgravity, there is extensive literature on 
mass and heat transport in oscillatory flows [26-36], the effect of gravity modulation on the 
classical B6nard problem [74-77] and the linear stability of planar fluid surfaces [78,79]. 

3.1 Order of magnitude analyses 

An order of magnitude (OM) or dimensional analysis of a physical process involves 
the examination of the system of equations that are assumed to govern the process. For 
fluid physics and material science experiments, the relevant physical processes are 
governed by equations describing the transport of heat, mass and momentum in single and 
multi-component fluid systems. The usual approach to OM analyses [51] involves the 
definition of reference quantities which appropriately characterize the physical system (i.e. 
reference scales for velocity, time, length, forces, etc.). Application of these scales to the 
(dimensional) governing equations, followed by the definition of characteristic 
dimensionless groups admits a comparison of the order of magnitude of every term in each 
equation. There are several comprehensive discussions of OM analyses applied to low- 
gravity situations (eg. [38,49-51]). 
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The advantage of the OM approach is that a great deal of information pertaining to the 
sensitivity of a given experiment can be obtained with little computational effort. The 
disadvantages are related to the fact that the choice of characteristic reference quantities for 
length, velocity, time etc. are not, in general, known a priori [80]. Thus, in cases for 
which the chosen characteristic scale is not (for a given set of conditions) representative of 
the actual process, the results of the analysis may be in error, often by several orders of 
magnitude. In addition, application of OM analyses to the problem of g-sensitivity has 
invariably involved implicit linearization and often poorly represents the multi-dimenional 
nature of the physical process. This can also lead to incorrect OM predictions. Finally, 
analyses are generally restricted to the examination of a single component disturbance (e.g. 
~ sin(cot) ). As a result, they may not properly indicate the response of the system to the 
typically complex multi-frequency disturbances characteristic of the spacecraft 
environment. It is nonetheless useful to examine the general trends predicted by such 
analyses since they can provide at least qualitative guidance for more detailed analytical and 
numerical studies. 

The determination of tolerable g-levels using OM analyses is based upon estimations 
of the response of a system to either steady accelerations or time-dependent disturbances of 
the form 

g(t) = a cos(cot), 

where co is the angular frequency of the disturbance. The susceptibility of an experimental 
system to such disturbances is defined via the magnitude of a particular response (say a 
temperature, compositional or velocity fluctuation) which must not be exceeded in order to 
bring the experiment to a successful conclusion. The most obvious trend predicted by 
analyses to date is shown in Figs. 3 - 5. The curves depict the maximum tolerable residual 
acceleration magnitude as a function of frequency, /= co/2jc, for given experiments and 


have the form 
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gmax — F(0),ncn t , Pi) . 

where licit is the magnitude of the maximum allowable response, and the pi, i=l,N 
represent the N material properties of the system (such as viscosity, thermal and species 
diffusivities). In genearl, below 10* 2 to 10*3 Hz the response is expected to be either 
independent or weakly dependent on frequency. At higher frequencies the analyses indicate 
that the allowable residual acceleration magnitude increases as approximately the square of 
the frequency. This behavior can be understood upon consideration of the following OM 
estimates [38] (used to obtain Fig. 3) for velocity, temperature and concentration changes 
associated with the acceleration (4): 

V _ g* A P 7 P (6) 

8 (co 2 + v 2 /L 4 ) 1 ' 2 


ST ~ V g 


A1ZL 

(co 2 + k 2 /L 4 ) 1/2 


(7) 


SC ~ V g 


AC/L 

(co 2 + D 2 /L 4 ) 1/2 ’ 


( 8 ) 


where g* is the residual acceleration magnitude, p is the mass density, and Ap/p is the 
change in density due to temperature and/or composition; AT/L, AC/L and L are, 
respectively, a characteristic temperature gradient, concentration gradient and length; k is 
the thermal diffusivity. The sensitivity limits are obtained by using either 8T, 8C or Vg as a 
sensitivity parameter and setting g* equal to the allowable residual acceleration gmax- 
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For situations which involve mixed thermocapillary convection and buoyancy-driven 
convection the governing dimensionless groups are [51]: the Bond number, B 0 = pg*L 2 /y, 
the Grashof number Gr = AT IPx I g*L 3 /v 2 > the Prandtl number Pr = v /k, and the Surface 
Reynolds number Rs = I y T IATL/jj.v. Here I y T I is the rate of change of surface tension y 

with temperature, px is the coefficient of thermal expansion, v is the kinematic viscosity, \x 
is the shear viscosity These groups respectively represent the ratios of buoyancy force to 
surface force, buoyancy force to inertial force, thermal and momentum timescales, and 
surface force to inertial force. 

The relative importance of gravity and theimocapillary forces can then be estimated 
upon considering the relative magnitudes of the velocity due to buoyancy and the velocity 
due to thermocapillary effects [50,51], The latter is given by 


V* = 


NAT 


(9) 


For thermocapillary flow to dominate, the ratio of Vg to V* must be less than one. 
For this to occur, the maximum magnitude of the acceleration must satisfy: 


g* < IN 


(co 2 + v 2 /L 4 ) 1/2 
vPt 


GO) 


Equation (10) can be re-interpreted in terms of the ratio of the Grashof number, Gr, 
and the surface Reynolds number Rs. The condition for buoyancy and thermocapillary 
forces to be the same magnitude is 


Gr = s(St 2 + l) 1 / 2 . 


( 11 ) 
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where St= cdL 2 /v is the Strouhal number. This condition is illustrated in Fig. 6. Three 
regimes of interest are identified. These regions are defined by the relative values of the 
Grashof, Strouhal and surface Reynolds numbers. In region I, the Strouhal number St is 
less than one, and the condition that thermocapillary forces dominate is equivalent to that 
for a steady flow. In region II, St is greater than one. In this region, for a fixed value of 
Rs, the value of Gr required to give buoyancy forces equal weight to thermocapillary forces 
increases with increasing values of the Strouhal number. This essentially reflects the fact 
that the characteristic time for the fluid response greatly exceeds the period of the 
disturbance. Thus, thermocapillary forces dominate at higher values than in region I. In 
region in buoyancy forces predominate over thermocapillary forces. For situations here a 
boundary layer scaling is appropriate i.e Rs » 1 [50]. The relative importance of 
buoyancy can be estimated via the ratio 


Gr = Rs 5/3 (St 2 + 1) 1/2 . (12) 

Figure 7 depicts the relationship between the estimated tolerable g-level and the 
frequency of the acceleration for experiments with the physical properties given in Table 2. 
These figures represent estimates based on two extreme flow regimes, penetrative or 
"viscous” flow (Fig. 7a) and boundary layer flow (Fig. 7b) d may be respectively 
considered as "worst" and "best" case estimates for the given sensitivity parameters. The 
curves have been determined for a sensitivity criterion which requires the tolerable 
acceleration levels to be such that the buoyancy forces should not exceed 10% of the 
thermocapillary forces (if the requirement had been 1%, the tolerable accelerations would 
be an order of magnitude smaller). In terms of acceleration levels predicted for the space 
station or measured on past missions the practical sensitive ranges are restricted to 
ldisturbances with frequencies less than 10 Hz. The majority of the estimates suggest that 
for these experimental parameters, periodic vibrations will generally not lead to significant 
buoyancy effects. 
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The above examples estimate the effects of simple single component accelerations. In 
reality, low-gravity disturbances tend to be associated with more than one frequency. 
Given that a system may respond to a multi-frequency disturbance in an "additive" way, 
tolerance curves such as these may underestimate the response to a given low-gravity 
disturbance. For example, models of the DMOS experiment, show that mass transport 
has an additive response to the (multi-frequency) residual acceleration [39]. Post flight 
analysis of the experiment results has demonstrated that the amount of mixing observed 
between organic liquids can be explained by the additive response of the system to a 
multicomponent disturbance. 

In many cases an a priori choice of length or time scales may not be obvious. For 
example, in a study of dopant (solute) uniformity in directionally solidified crystals, the 
OM estimates of Rouzaud et al. [52] and Camel and Favier [83] are in agreement with the 
direct numerical simulations of Chang and Brown [84] for a Schmidt number (ratio of the 
melt's kinematic viscosity to dopant diffusivity) of fifty but overestimate the amount of 
radial segregation for a Schmidt number of ten. 

The results of an order of magnitude or scaling analysis have been compared in [80] 
with those of numerical simulations of the effects of steady residual acceleration on 
compositional non-uniformity in directionally solidified crystals. The basic model consisted 
of a Bridgman-type system with an ampoule translated between hot and cold zones 
separated by a thermal barrier. In order to apply the Camel-Favier technique [52,83] to the 
numerical simualtion, it was necessary to evaluate a proportionality constant, x 0 2 . This 
relates the maximum convective velocity obtained in the numerical simulations with the 
Grashof number. The value of this constant depended on the orientation of the acceleration, 
and on the chosen length scale. The adiabatic zone length L a yielded results which best 
showed the trends predicted by the Camel-Favier approach. Table 3 lists values of x 0 
associated with different choices of L a as a reference length. Having calculated Xq we then 
determined the transport regime (see references [52,80,83] for details) by graphing, in Fig. 
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8, GrSc vs Peg for the case of the residual acceleration parallel to the interface. Here Sc 
(Schmidt number) gives the ratio of kinematic viscosity to thermal diffusivity, and Pe g = 
VL /D, is a ratio of the translation rate to the characteristic diffusive speed where L is the 
characteristic diffusion length and D the solute diffusivity. The lateral composition non- 
uniformities in the crystal obtained from numerical simulations [9] are also given in Fig. 8. 
A comparison between computed and predicted non- uniformity values showed that at high 
values of GrSc (region Ic in Fig. 8), the order of magnitude approach meets with some 
success but is less faithful at low values of GrSc (for example the transition between 
regions lb and II). Fig. 8 also contains two cases for which Pe g is reduced at fixed GrSc. 
According to the estimates the lower GrSc case should have yielded a higher compositional 
non-uniformity as the growth rate was reduced, i.e. the estimates predict that the system 
moves closer to the convective diffusive transition. This was not reflected in the non- 
uniformity obtained from the numercial simulations. Evaluation of the numerical results 
showed that the increase in the characteristic diffusion length scale (associated with 
reduction in growth rate) led to an order of magnitude decrease in the solute gradient at the 
interface. Hence, because the system-wide variation in composition is small, the interfacial 
non-uniformity is correspondingly small. 

The above results indicate that scaling arguments alone can at best be expected to 
yield order of magnitude accuracy. The limitations of scaling arguments are hardly 
surprising if one considers the neglect of the multi-dimensionality (boundary conditions, 
flow structure, etc.) of the physical situation which is inherent in such scalar descriptions. 
Of course, prior knowledge of the system behavior can be used to locally improve the 
accuracy of the scaling, (eg. Kimura and Bejan [85]). Similarly, if a system locally exhibits 
one-dimensional behavior, such as boundary layer flow, then the appropriate scaling for 
the formal reduction of the transport equations can be used effectively to estimate the local 
system response. In general, however, our understanding of the response of systems to 
residual acceleration can be furthered only by an approach which uses scaling. 
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mathematical modelling and, naturally, the results of experiments in a complementary 
fashion. 

3.2 Analytical models 

A few attempts have been made to assess the effects of time-dependent residual 
acceleration using analytical models [53-59]. The first attempt to model the effects of 
disturbances on heat transfer between a fluid and a solid bounding medium was carried out 
by Gebhard [53]. Under the assumption that the fluid and its container would be subject to 
a "sequence of abrupt relative displacements spaced by a time interval x c , •••". Gebhard 
took T c to be a random distributed variable. He showed that as T m (the most likely value of 
x c ) decreases, the resulting heat transfer was greater than predicted by conduction in the 
absence of a disturbance. In particular, it was found that for T m < d 2 /ic (where d is a 
characteristic length and k is the thermal diffusivity) the heat transfer increased rapidly as 
Tm decreased. The main limitation of this model is that fluid convection (and thus fluid 
properties such as kinematic viscosity) do not enter the analysis. It is well known that the 
response of heat transfer to convection varies according to the Prandtl number. In 
particualr, for low-gravity flows (as we shall see later) the temperature field is less sensitive 
to convection for low Prandtl number fluids (Pr < 1 ). 

Approximate solutions for transient convection in a cylinder with an azimuthal 
temperature variation have been obtained by Dressier [54], Residual accelerations 
representing the motion of an astronaut and a transient rotation of the spacecraft were 
imposed on the system. Linear accelerations were imposed perpendicular to the cylinder 
axis. For kinematic viscosities of 10* 2 cm 2 s**, cylinders of radius 1 and 2 cm with a 
maximum temperature difference of 95 K across the diameter were examined. The 
maximum fluid velocities consequent to the astronaut motion (modelled by a sequence of 
two step functions separated by 2 seconds with 4 x 1CH g magnitude, 0.5 second duration 
and zero mean, oriented perpend’ ~ular to the temperature gradient) were 12 (im s' 1 for 
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both the large and small diameter cylinders (see Fig. 9). For the larger cylinder the second 
pulse reduced the velocity from 12 pm s' 1 to 0 by the end of the pulse. The effect of the 
second pulsewas to change the velocity in the small container from 12 to -3 pm s' 1 by the 
end of the pulse. The velocity decayed asymptotically to zero reaching a value close to zero 
after 15 seconds. The reason the response was different for the two containers can be 
explained simply in terms of the momentum decay times (ti/2 ~ R^/v) for each container. 
The velocity in the larger container had hardly begun to decay when the second impulse 
arrived. The shorter decay time associated with the smaller container meant that the unpulse 
acted on a smaller magnitude velocity which explains the change in direction of the motion. 
While the disturbance had a mean value of zero the fluid particles attained a finite 
displacement (indicated by the value x of the angular displacement in Fig. 6). Similarly, a 
disturbance corresponding to a simple spacecraft maneuver was examined. The simulated 
maneuver corresponded to a 90° rotation of the spacecraft at a rate of 1° s The 
maximum velocities associated with the centrifugal acceleration were, respectively, 40 and 
35 pm s' 1 for the large and small diameter cylinders. 

Experiments involving free surfaces of liquids have a high probability of responding 
to residual accelerations. Indeed, unexpected results have been attributed to residual 
accelerations in several cases [40,86-88]. Particular examples of such experiments involve 
liquid (float) zones. Three types of liquid zones are typically the subject of low-gravity 
experiments [89]: isothermal and non-isothermal zones suspended between inert solids (i.e. 
liquid bridges), and non-isothermal zones suspended between a feed rod and a growing 
crystal. The dimensionless group associated with the shape of isothermal static liquid zones 
is the Bond number, B 0 , defined earlier. It expresses the relative importance of 
gravitational and surface forces. A comprehensive description of stable zone shapes in zero 
gravity, and as a function of B 0 , is given by Martinez et al. [41]. Fig. 10 displays then- 
results for the stability of static zones with a steady acceleration parallel to the zone axis. In 
order to relate the Bond number to a given acceleration level they give the value of B 0 = 1.4 
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for a 1 liter volume of water subject to an acceleration of 1 cm s' 2 . For B 0 = 0, the zone 
volume V=V*= rcR 2 L and contact angle a = 90° the maximum stable length is 2 tiR (the 
Rayleigh limit). Coriell et al. [90] investigated the effect of B 0 on the maximum length of 
such zones and compared theoretical predictions with experimental results (see Fig. 11). 
The actual length of the equilibrium zone is determined by the volume, the contact angle a 
and the static Bond number. 

Estimates of the effect on the free surface shape of oscillatory residual acceleration 
parallel to the axis of an isothermal liquid bridge have been made by Langbein [55]. Not 
suprisingly, his calculations suggest that the sensitivity of the surface is highest for 
disturbances with frequencies close to the natural frequency of the bridge, typically in the 
0.001-10 Hz range. For the example given in [55], the maximum acceleration magnitudes 
that can be tolerated by a bridge less than 90% of the maximum stable bridge length (given 
by the Rayleigh limit) range from 5xl0' 8 to 5x10' 3 g. The criterion used to determine the 
tolerable residual acceleration was that the amplitudes of all surface shape excitations with 
frequencies other than the resonance frequency be kept to 0.1% of the radius of the 
column. Examples of the results of this analysis are shown in Fig. 12. The curves 
represent the maximum tolerable axial acceleration for a given aspect ratio (Fig. 12a), and 
for a surface response with a given number of nodes (Fig. 12b). The results express the 
maximum tolerable frequency ao in terms of the g-jitter frequency to, the resonant 
frequency, co m , of the liquid zone associated with the m ( h node of the deformation of the 
zone surface, the viscosity of the liquid v and the spatial wavenumber 

ao < min[V (co 2 -co£) - covq m ] , (9) 

where c4 = yq^[(q m R) 2 - l]/pR , and q m = (m+l)tt/H, R is the column radius and H 
the column height. Note that the global minimum, and local minima of the tolerable 
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acceleration curve occur at disturbance frequencies equal to the natural frequencies of the 
system. 

Zhang and Alexander [91] have also analyzed this problem using a slender-body 
approximation. The problem of determining the axisymmetric response of the shape of the 
free surface of a cylindrical liquid column bounded by two solid regions is modelled by a 
1-D system of non-linear equations. It is found that the sensitivity of the zone to breakage 
and shape changes depends on the static Bond number, aspect ratio and viscosity as well as 
the amplitude and frequency of the disturbance. The general trend is an increase in 
tolerable residual acceleration with increasing frequency. At the eigenfrequencies of the 
zone, however, there are strong deviations from this trend. At these frequencies this model 
also shows that the tolerable residual gravity level can be two orders of magnitude lower. 
These most sensitive frequencies have been found in the 10*^- 10'^ Hz range (see for 
example Fig. 13). For the cases examined, maximum tolerable residual gravity levels as 
low as 10^ g have been calculated. For higher viscosities the tolerable acceleration level is 
increased for all frequencies. The equilibrium shape, as determined by the steady 
background acceleration, has a prounounced effect on the response at low frequencies. A 
change in slenderness of the bridge markedly changes the sensitivity to residual acceleration 
as the Rayleigh limit is approached. It is clear that for liquid bridges, small magnitude 
accelerations cannot be ignored. 

For non-isothermal zones, the dimensionless groups introduced earlier in section 3.1 
govern the problem. For molten silicon Sekerka and Coriell [89] have calculated the values 
of these groups. They are presented in Table 4. The relative magnitudes of the terms 
suggest that on earth the zone shape will be dominated by the capillary effect (large B 0 ) for 
typical zone sizes, while under low-gravity conditions the zone shape will be influenced 
mainly by the temperature gradient along the interface. In addition (as discussed earlier) 
dynamic distortions may also be important. Considerations of the float zone crystal growth 
process introduce further complications which ultimately require independent analysis. In 
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the previous cases, the surfaces of the inert solids between which the zone is suspended are 
usually planar. This is generally not the case for the crystal growth process. The actual 
shape depends primarily on the thermal conditions. The influence of the meniscus angle at 
either of the crystal-melt-vapor tri-junctions can also be important and has been analyzed by 
Heywang [93] and Coriell et al. [94]. In addition to zone failure due to the classical 
capillary instability (or Rayleigh instablity), a "dewetting” instability is found to be 
important. The instability arises when the angle between the meniscus and the melt-crystal 
interface is less than the wetting angle for the liquid melt on the solid crystal. For Si zones 
of greater than 1 cm radius, zone failure on earth occurs due to the Heywang instability 
[94,95]. Zone failure by this mechanism is unlikely under low-gravity conditions. 

Feuerbacher et al. [56] developed an ad hoc model to describe experiment sensitivity 
to residual accelerations. The model is based on the motion of a body in response to an 
applied sinusoidal oscillation. The governing equation is that of a damped harmonic 
oscillator. The tolerable residual acceleration are (as expected) found to increase with 
increasing frequency of the disturbance. Furthermore, there is a minimum tolerable 
acceleration at disturbance frequencies corresponding to the natural frequency of the 
system. 

The effect of time dependent accelerations on experiments involving drops and 
bubbles has been examined for the cases of sinusoidal and Heaviside-function disturbances 
[57]. The criterion used to define the sensitivity is that the displacement of the drop not 
exceed 10% of the drop radius. For a given set of physical properties, the admissible 
accelerations can be expressed as a function of frequency (Fig. 14). The tolerable 
acceleration levels exhibit trends similar to those depicted in Figs. 4-6. 

The influence of the gravity gradient, atmospheric drag and spacecraft attitude motion 
on the motion of a spherical solid particle immersed in viscous fluid has also been 
considered [5,6]. The extent and nature of the influence of the gravity gradient and the 
atmospheric drag on the trajectory of the particle are, for a given value of the Stokes drag 
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coefficient, shown to be dependent on the characteristic orbital attitude, the magnitude of 
the atmospheric drag accelerations and the distance of the particle from the mass center of 
the spacecraft. 

Amin [59] has investigated the influence of g-jitter on heat transfer from an isothermal 
sphere maintained at a constant temperature greater than the ambient fluid temperature. The 
body force was taken to have a single frequency sinusoidal time dependence and a small 
amplitude. A significant result of the analysis is that the Reynolds stresses associated with 
the fluctuating flow result in a steady streaming motion. This has implications for both heat 
and mass transport in oscillating flows. For the low Pr fluids examined, it was concluded 
that buoyancy induced convection caused by high frequency g-jitter (with amplitudes 10‘ 2 
g) will result in heat transfer conditions significantly different from pure conduction. 
However, it was found that in high Prandtl number fluids with small kinematic viscosities, 
low frequency g-jitter will influence on heat transfer. 

Analyses of heat transfer in (laminar) oscillating flows in cylinders and between 
parallel plates [27-33] have shown that at high frequencies the effective thermal diffusivity 
is proportional to 

KelT~ t!i ^£ (10) 

W 0 

while for low frequencies 


Keff« co Ax 2 W 0 . 


(ID 


Here co is the circular frequency. Ax is the cross-stream average displacement of a fluid 
element over half the period of the oscillation and W Q = L(oVv) 1 ^ 2 is the Wormsley number 
which represents the ratio of the viscous response time to the period of the oscillation. For 
a given frequency the heat transfer process was found to be most efficient for Pr = Jt /W 0 2 
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i.e. when the characteristic time, = L 2 /k associated with heat diffusion is equal to half 

the period of the oscillation. Fig. 15 illustrates this result for fluids with different 
diffusivities and L = 1 cm. 

3.3 Linear stability and gravity modulations 

The stability of a heated fluid layer of infinite lateral extent is affected by gravity 
modulation and has received a limited amount of attention for the case of sinusoidally 
modulated gravity [74-77]. The static case, where the gravity vector is perpendicular to the 
layer and parallel to the temperature gradient, is stable unless the critical Rayleigh number is 
exceeded. The effect of gravity modulation is to introduce stability for given combinations 
of modulation amplitude and frequency for situations which in the absence of modulations 
would exhibit instability. In contrast, for the case where the gravity vector and the 
temperature gradient are anti-parallel (which is stable for a steady gravity vector,) 
instability can occur due to gravity modulation. 

The work of Biringen and Peltier [77] focused on the behavior of a 3-D Rayleigh 
Bdnard system subject to spatial and temporal gravity modulation. They used a pseudo- 
spectral method to examine the system numerically and studied the effects of sinusoidal and 
random gravity modulations with the gravity vector oriented perpendicular to the fluid 
layer. It was found that when the base acceleration was reduced from lg to zero, the 
system parameters that had resulted in a synchronous response at 1 g led to a conductive 
state for the modulated zero g case. The spatially random modulations about a zero g base 
acceleration led to excitations of the 3-D temperature field. These excitations were less 
pronounced than for the 1 g case. For temporally random modulations about zero g, the 
cases considered that where excited at 1 g no longer exhibited excitation. 

The effect of oscillating accelerations and impulse forcing of an interface separating 
two immiscible fluids has recently been considered [78,79]. Both these analyses are 
restricted to time-dependent accelerations with zero mean which are oriented perpendicular 
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to the interface between the fluids. The linear stability of a basic state characterized by a 
planar interface and an oscillating pressure. For the case of an arbitrary oscillating 
disturbance it is shown that the linear evolution of an infinitesimal perturbation of the 
interface is governed by a single integro-differential equation. Except in the limit of zero 
surface tension the effect of viscosity is shown to be small. As surface tension becomes 
dominant, resonant instabilities are excited at lower and lower wavenumbers. The interface 
is found to be more unstable to impulse type forcing than to sinusoidal forcing [79]. 

3.4 Direct numerical simulation 

The effects of low-gravity on the transport of heat and momentum have been 
examined in a number of articles [37, 60-71]. Robertson et al. [60,61] found that for 
convection in circular cylinders with azimuthal variations in temperature at the boundary 
and the gravity vector applied perpendicular to the cylinder axis, the intensity of convection 
follows the prediction of Weinbaum's first order theory [96] for low Rayleigh numbers. 
Weinbaum's first order theory predicts a simple sinusoidal dependence of the maximum 
velocity as a function of orientation of the gravity vector. The effects of a variety of 
acceleration vectors (sinusoidal, cycloidal and linear periodic) on motion in three fluids 
corresponding to mercury, helium and water (10" 2 < Pr < 10) were studied by Spradley et 
al. [62], They found that the steady mean part of the applied disturbance is more important 
than the oscillatory part (frequency =1 Hz) in determining the flow field and heat transfer. 
Kamotani et al. [37] solved a linearized approximation of the Boussinesq equations and 
investigated the effect of an applied acceleration consisting of a time mean part and an 
oscillatory part on the temperature and flow fields in a rectangular enclosure. They found 
that the thermal convection was predominantly oscillatory in nature. In addition they found 
that the acceleration perpendicular to the temperature gradient is the most important for the 
generation of fluid motion. 

The response of natural convection of a Boussinesq fluid in a cylinder to residual 
acceleration has been examined by Heiss et al. [64], and Schneider and Straub [97], In the 
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work of Heiss et al. the calculations are three-dimensional and involve gravity pulses and a 
rotating gravity vector. The ends of the cylinder are held at different temperatures. The 
cylinder walls are taken to be either perfectly insulating or perfectly conducting. The results 
are presented in terms of the relevant dimensionless parameters in the original papers. 
Figure 16 illustrates their results for impulse disturbances. The physical conditions were 
taken to correspond to: a cylinder of aspect ratio 1, with diameter = 1 cm, v = 5x10' 3 cm 2 
s' 1 , K= 7. 14x10" 3 cm 2 s' 1 , (3 = 2.5xl0' 4 , and a temperature difference, AT, of 25° K. 
The magnitude of the acceleration pulse (duration =0.7 seconds) corresponding to a 
Rayleigh number, Ra*, of 2000 is ~ 3x10*3 g results in a maximum velocity of 
approximately 5.8 cm s' 1 . The decay time for the flow field is about 5.6 seconds. Notice 
that for all the results, the decay time of the pulse is independent of the magnitude of Ra*. 
For a rotating gravity vector of a given magnitude, the maximum velocity of the system is 
independent of the dimensionless frequency / = / *d 2 /k, up to / ~ 10 whereupon it 

decreases rapidly with increasing frequency. For the fluid with the physical properties 
described in the previous example, / * corresponds to 7.14xl0' 2 s' 1 . 

Schneider and Straub [97] examined the effect of Prandtl number on the response to 
g-jitter for the same system as Heiss et al. They examined the response of the system to 
pulses of infinite length and sinusoidally varying pulses. For sinusoidal pulses they found 
for systems with Pr =0.71 and 7.0 (air and water) that the maximum velocity decreased 
with increasing frequency for dimensionless frequencies, F=/L 2 /k, greater than 10. At 

Pr=7 a maximum in the velocity was observed between F = 10 and 100 for the maxium 
amplitude pulse examined (Ra = 2xl0 5 , or I g I = 3x1 O' 4 g for a fluid with the properties 

of water, a 1 cm diameter cylinder and a temperature difference of 10 K). They determined 
that the system was least sensitive when the transient accelerations acted along the cylinder 
axis. 

The extent to which gravity causes buoyancy-driven fluid motion (and thus, solute 
redistribution) during directional solidification has been examined using numerical models 
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of buoyancy-driven convection in cylindrical and rectangular geometries[9,65-69, 84,98- 
101]. Except for [9,65-69] these studies are restricted to axisymmetric situations in which a 
steady gravity vector is oriented parallel to the axis of a cylindrical ampoule. McFadden 
and Coriell [66] have undertaken 2-D calculations of the effects of time-dependent 
accelerations on lateral compositional variations during directional solidification. The 
gravitational acceleration was assumed to have a uniform magnitude and rotation rate. The 
amount of compositional non-uniformity was found to increase with decreasing rate of 
rotation (see Fig. 17). 

Two- and three- dimensional models of directional solidification from dilute gallium- 
doped germanium melts have been used to determine the sensitivity of crystal 
compositional uniformity to both time-dependent and steady residual accelerations 
[9,67,69]. The specific boundary conditions, thermo-physical properties of the melt, 
growth rates and ampoule size are all found to play a role in the determination of the 
experiment sensitivity. For a given set of operating conditions, it is found that at growth 
rates on the order of 6.5 mm s' 1 the orientation of the experiment with respect to the steady 
component of the residual gravity is a crucial factor in determining the suitability of the 
spacecraft as a means to suppress or eliminate unwanted effects caused by buoyant fluid 
motion. The worst case appears to be when the acceleration vector is parallel to the crystal 
interface. At growth rates on the order of microns per second, this orientation leads to 
compositional non-uniformities of 10-20% when the magnitude of the acceleration is of the 
order 10‘ 6 g. If, however, the growth rate is lowered by an order of magnitude, the non- 
uniformity is reduced significantly (down to 1- 5% in this case). Table 4 summarizes the 
results obtained for steady accelerations. 

The directional solidification (Bridgman-Stockbarger) process is also extremely 
sensitive to transient disturbances. The response of the system to a variety of impulses will 
be discussed. For example, a 3x1 0"3 g impulse of one second duration acting parallel to the 
interface of a growing crystal produces a response in the solute field which lasts for nearly 
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2000 seconds. Consequently lateral and longitudinal compositional variations occur over a 
length of nearly 6 mm in the grown crystal. 

The response of the solute field and the lateral non-uniformity to oscillatory 
accelerations varies from no response at all (at frequencies above 1 Hz with amplitudes 
below 10* 3 g) to a significant response at 10' 3 Hz at amplitudes on the order of 10* 6 g. In 
addition, additive effects were observed for combinations of a steady and a low frequency 
residual acceleration component. These additive effects gave rise to significant lateral and 
longitudinal non-uniformities in concentration. A number of different types of periodic 
disturbances were employed. Single frequency disturbances of the form g(t) = go + gn 
cos(2rc/ n t) were examined with g Q = 0, V2(10)* 6 and V2(10)' 5 g, oriented parallel, 
perpendicular and at 45° to the crystal-melt interface. The frequency range examined was 
f n = 10 * 4 , 10' 3 , 10* 2 , 10* 1 , 1 and 10 Hz. For frequencies greater than 10' 2 Hz, there were 
no discemable effects on the solute fields. The velocity field did, however, respond to the 
oscillatory disturbances. For the case of 10" 3 Hz (at 5xl0* 6 g) the response of the solute 
field was significant. Lateral and longitudinal non- uniformity levels in excess of 15% were 
calculated. Figures 18-20 show the lateral non-uniformity as a function of time and 
highlights the additive effect of oscillatory and steady components of the residual 
acceleration. 

The effect of a multiple frequency disturbance is illustrated in Figure 20. The 
acceleration consists of steady and periodic contributions with the form: g(t) = go + gi 
cos(2ti10' 3 t) + g 2 cos(2tt10' 2 t ). Here I g 0 1 = V2xl0* 6 g, I gi I = 3V2xlO* 6 g and I g 2 I 
= 3^2x1 0' 3 g. The magnitude of the compositional non-uniformity £ is seen to vary with 
the frequencies of the acceleration. More recently, it has been shown [67] that, particularly 
at high frequencies, it is important to consider the long time behavior of the system, 
because for high frequency disturbances ( * 0.1 - 1 Hz) the transient behavior of the 
system is more sensitive than the long time behavior. For example, the response to a 1 Hz 
10' 2 g acceleration was found to exhibit a concentration non-uniformity on the order of 
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10% during the first 80 seconds whereupon it decayed slowly to 1% of over a period of 
2500 seconds. In addition, calculations have been made using samle acceleration time 
histories obtained on Spacelab 3 (SL-3). The results varied according to the type of SL-3 
data that were input. For the example shown in Fig. 21. the response of the solute field 
was negligible, although the velocity field fluctuated on the order of the crystal growth rate. 
The system exhibited some response to transient disturbances (for example the structural 
response to thruster firings) that were superimposed on background low frequency 
acceleration. Low frequency (< 10' 1 Hz) acceleration components with amplitudes of 10' 3 
g or more were also found to give rise to undesirable non-uniformity. 

Griffin and Mokatef [72,73] modelled the nature of melt convection in a Bridgman 
system subject to steady and unsteady axial acceleration. They also found that at conditions 
corresponding to low-gravity, the response time of the melt velocity is proportional to the 
momentum diffusion time (i.e. it is controlled by the characteristic system length scale and 
melt viscosity). In accord with the general trend exhibited by most physical systems 
discussed in this article, they found that the velocity response to sinusoidal g-jitter 
decreased as the inverse square of the momentum diffusion time. 

A general conclusion that can be drawn from all attempts to characterize gravity- 
driven convective effects on directional solidification from two component melts is that the 
maximum lateral solute non-uniformity (radial segregation for the axisymmetric cases) 
occurs near the transition from diffusion dominated to convection dominated growth 
conditions [52,84,9]; that is, when convective velocities are of the same order of magnitude 
as the diffusive velocities. The conditions under which this "transition" takes place will 
depend on the specific nature of the forces driving convection. The orientation of the steady 
component of the gravity vector is crucial in determining the magnitude of the gravity 
vector at which this transition occurs. Thus, for a given set of operating conditions the 
orientation of the gravity vector determines the suitability of a low-gravity environment for 
directional solidification experiments. 
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Motivated by the growth of triglycine sulfate (TGS) crystals on Spacelab 3, [11], 
Nadarajah et al. [102] have examined the acceleration sensitivity of crystal growth from 
solution. Thermal and solutal convection were included in their analyses and the crystal 
growth rate was chosen as the sensitivity parameter for the response to convective 
transport. Simulations were carried out for steady, impulsive and periodic accelerations in 
order to determine tolerable acceleration levels. Long-time simulations of the experiment 
were conducted with steady background accelerations of 10~ 6 and 10‘ 5 g. Impulsive and 
periodic disturbances of higher magnitudes were imposed at intermediate points. For 
steady accelerations the system was shown to have diffusion dominated transport at 10' 6 g 
but for 10' 5 g a transition to convection dominated transport occurs. Figure 22 summarizes 
the results for oscillatory disturbances. For each case, the numbers correspond to the ratio 
g|3A c/f v*, which is an estimate of the the ratio of buoyancy to inertial forces. Here [5 is the 
solute expansion coefficent, Ac the characteristic concentration difference, / is the 
frequency of the acceleration and v* is the calculated characteristic fluid velocity. Clearly, 
for this example, a simple estimate alone will not suffice to determine the system 
sensitivity. The system is relatively stable to impulsive and periodic disturbances unless 
their magnitudes are very large. 

4. Summary and discussion 

The work described in this paper examines the sensitivity of a variety of space 
experiments to residual accelerations. In all the cases discussed, the sensitivity is related to 
the dynamic response of a fluid. In some cases the sensitivity can be defined by the 
magnitude of the response of the velocity field. This response may involve motion of the 
fluid associated with internal density gradients, or the motion of a free liquid surface. For 
fluids with internal density gradients, the type of acceleration to which the experiment is 
sensitive will depend on whether buoyancy-driven convection must be small in comparison 
to other types of fluid motion (such as thermocapillary flow), or fluid motion must be 
suppressed or eliminated (such as in diffusion studies, or directional solidification 
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experiments). In the latter case, the experiments (for example diffusion experiments and 
directional solidification experiments) are sensitive to steady and low frequency 
accelerations (< 10'^ Hz). For experiments such as the directional solidification of melts 
with two or more components, determination of the velocity response alone is insufficient 
to asses the sensitivity. The effect of the velocity on the composition and temperature field 
must be considered, particularly in the vicinity of the melt-crystal interface. As far as the 
response to transient disturbances is concerned the sensitivity is determined by both the 
magnitude and frequency the acceleration and the characteristic momentum and solute 
diffusion times. 

The directional solidification of doped germanium crystals has been carried out on 
Skylab [17,18,103,104] and the Apollo-Soyuz Test Project (ASTP) [105,106]. The latter 
involved the directional solidification of seeded Ge melts. Melts were doped with gallium 
and some were doped with 1% Si and .001% Sb [105]. The results of these experiments 
revealed strong asymmetric non-uniformities in the space grown crystals. Lateral variations 
were also observed in samples grown under terrestrial conditions but were much less 
pronounced [25]. It has been argued that the asymmetric redistribution of the dopant can be 
ascribed to "barometric diffusion" of the solutes due to the acceleration gradient in the melt 
arising from the rotational motion of the spacecraft [105,106]. The basis of the argument, 
however, appears to ignore the presence of gravity gradient and atmospheric drag effects 
and does not explicitly account for the spacecraft attitude motions. It has been demonstrated 
[9] that melt convection can occur in response to a low magnitude steady residual 
acceleration, and that whenever the acceleration vector is not aligned with the ampoule axis 
strong asymmetries in composition can occur. Since there is a non-linear dependence of 
compositional uniformity on convection, there is always the possibility that the observed 
non-uniformities in the terrestrially grown samples will be smaller than those observed in 
space grown samples. Thus, convection caused by residual acceleration cannot be 
discounted as the origin of the non-uniformities. 
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Thermocapillary flow experiments require that buoyancy effects are small in 
comparison to capillary effects. The estimated range of accelerations to which they are 
sensitive lies between 5xl0'3 Hz and 2 Hz for the most sensitive cases. Unless the 
magnitude of the effective buoyancy force (expressed through the Grashof number) 
associated with a given acceleration is on the order of or greater than Rs(l+ St 2 ) 1 / 2 , 
buoyancy effects are not likely to be significant. It is possible that more subtle effects of 
the acceleration environment might occur (for example influencing the transition from 
steady to oscillatory flow) but this cannot be determined by order of magnitude analysis. It 
is more likely that thermocapillary flow experiments will be most sensitive to the motion of 
the free surface in response to g-jitter. 

The response of free surfaces to g-jitter has been examined through analytical and 
numerical models. For cylindrical liquid bridges this has been undertaken only for axial 
accelerations. The shape of the surface is most responsive to disturbances close to the 
natural frequency of the bridge and its harmonics. The sensitive frequencies lie in the 
5xl0‘ 3 10 Hz range, with maximum tolerable magnitudes from 10* 8 - 10* 3 g depending on 

the sensitivity criterion and the particular model employed. 

Figure 23 depicts the range of measured acceleration amplitudes (courtesy of Dr. H. 
Hamacher, DLR, and M. J. B. Rogers, University of Alabama in Huntsville) as a function 
of frequency for several “shuttle” missions. It should be remembered that in addition to 
the accelerations shown here there are gravity gradient and atmospheric drag accelerations 
with magnitudes on the order of 10*6- 10* 2 g which can be steady (for gravity gradient 
stabilized type attitudes) or have frequencies that are twice the orbital frequencies (for a 
solar inertial attitude). In Fig. 24 we depict these together with g-tolerance curves selected 
from the examples discussed earlier. There is some overlap between the sensitivity curves 
and measured accelerations. Even if the fact that most of the curves are obtained for order 
of magnitude estimates is taken into account, it is clear that a careful evaluation of the 
interaction between experiments and the residual acceleration environment is necessary. 
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Once the experiment sensitivity has been determined (which may sometimes require a more 
detailed analysis than any discussed in this review) it should be possible for the 
experimenter to choose operating conditions which will minimize the effect of the predicted 
acceleration environment. Since flight opportunities are limited, it is particularly important 
that the chances of success for any given experiment be greatly improved by detailed 
modelling prior to flight. This will allow investigators to better anticipate undesirable 
effects of the prevailing residual acceleration conditions, and will ultimately lead to a better 
overall understanding of the physical system. 
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Table 1. 


Changes in residual acceleration components [p.g m'^] along the local 
vertical (Xi) and perpendicular to the oribital plane (X 2 ) associated with the 
gravity gradient stabilized attitude, and estimated atmospheric drag 
accelerations, a and b are, respectively, best and worst drag estimates for 
the “Shuttle” (1) and NASA’s Space Station (2). From [5] 




Table 2. Physical properties and dimensionless groups for various experimental systems: 1- di methyl silicone oil [44], 2- methanol 
[81], 3- NaN 03 [45], 4 - indium [82], 5- selenium - [82], 6- silicon [15]. 
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Table 3. Comparison of numerical (£* m ) and order of magnitude estimates ( m ) of compositional 
non-uniformity for different choices of length scale [80], The non uniformity x in these 

cases is given by ^ | c(0) - c^ 1/ c(0), c^ is the melt composition far from the crystal 

and, because k< 1 , c(0) is the maximum value of the composition at the interface. 
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Table 4. Values of dimensionless groups for silicon calculated from the data of Chang and 

Wilcox [92], The length scale is taken to be the zone radius R; p T is estimated at 1.6 x 10- 

4 K 1 . For space, g s is taken as lO^g. AT is assumed to be 10 C*. (After Sekerka and 
Coriell [89]). 
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Table 5. 

Compositional 

non-uniformity % [%] for Ge:Ga [9,80]; £, = (c sma x - 

Csmin )/ c sav. 


where c s denotes the crystal composition at the melt-crystal interface. 
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e g is the unit vector parallel to g, N is the normal vector to the crystal. L = 1 cm for all cases. Values in 
parentheses indicate 3-D results. All calcualtions were undertaken using the Thermo-physical properties of gallium- 
doped germanium (9) and A) and B) refer to the operating conditions listed below. 
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Figure Captions 


Fig 1. 

Fig 2. 

Fig. 3 

Fig. 4 

Fig. 5 
Fig. 6 

Fig. 7 


An example of a relatively “quiet period” measured on SL-3 (< 4xl0' 3 g, g 
= 9.8 m s' 2 ), a) Total acceleration array, b) Combined amplitude spectrum, 
note components in the 5-6 Hz and 7-8 Hz, 10-12 Hz, 17Hz and 30-35 Hz 
ranges. After Rogers and Alexander [47]. 

“Window” containing three thruster pulses, a) Total acceleration array, 
maximum magnitude ~lxl0' 2 g. b) Combined amplitude spectrum, 
dominant frequency components 11 Hz and 17 Hz. After Rogers and 
Alexander [47]. 

Estimated tolerable residual accelerations (in units of g = 9.8 m s' 2 ) as a 
function of frequency for: 1) a fluid physics experiment involving a 
temperature gradient, 2) a crystal growth experiment and 3) a thermo- 
diffusion experiment. After Monti et al. [38]. 

Estimated tolerable residual accelerations (in units of g = 9.8 m s' 2 ) for 
semiconductor and metal solidification experiments. The sensitivity 
parameter is longitudinal segregation, Ac/c . Pe s = 5, temperature gradient = 
50K cm -1 , k=0.1. After Monti et al. [36], 

Estimated tolerable residual acceleration (in units of g = 9.8 m s' 2 ) ranges 
for diffusion experiments. After Monti et al. [36]. 

Tolerable acceleration limits for a thermocapillary flow experiment 
expressed in terms of Grashof (Gr), Strouhal (St) and Surface Reynolds 
(Rs) numbers. 

Tolerable acceleration limits (in units of g = 9.8 m s' 2 ) as a function of 
frequency for the thermo-capillary experiments with properties listed in 
Table 2. a) Gravity forces 1% of thermo-capillary forces, b) gravity forces 
10% of thermo-capillary forces. 


Fig. 8 


Fig. 9 


Fig. 10 


Fig. 11 


Fig. 12 


A plot of GrSc vs. Peg for the results of numerical simulations for steady 
accelerations parallel to the crystal-melt interface [80]. The solid and dashed 
lines delineate the following transport regimes: I a - c convective regime; II 
advective-diffusive regime [83]. Note that the adiabatic zone length has been 
used to calculate the Grashof and growth Peclet numbers, and that x 0 
=0.27. Numbers in parentheses indicate the value of compositional non- 
uniformity £ = lc smax - c sm j n l/c sa v x 100%, where c s denotes the crystal 
composition at the melt-crystal interface, associated with each of the 
numerical results. 

Transient convection in large and small test cells due to residual 
accelerations caused by astronaut motion after Dressier [54]. Here x is the 
angular displacement, v the velocity and ti/2 is the decay time given by 
0.68xR2/v. 

Stability limits of an isothermal liquid bridge of volume V and length L 
suspended between discs of radius R and subject to a steady axial 
acceleration defined by a Bond number B 0 , after Martinez et al. [41]. The 
angle 0 is defined in the inset 

The maximum stable zone length of an isothermal liquid zone suspended 
between inert solids of radius R as a fucntion of Bond number: Comparison 
of theory and experiment. After Coriell et al. [90]. 

Tolerable residual accleration (in units of g = 9.8 m s' 2 ) for experiments 
investigating resonance frequencies of liquid columns subject to axial 
acceleration , following Langbein [55]. a) sensitivity for different aspect 
ratios A = L/2jtR, b) change in sensitivity with number of nodes in the 
deformation for A = 0.9999. 



Fig. 13 


Curves of tolerable acceleration (in units of g = 9.8 m s’ 2 ) versus frequency 
for a maximum liquid bridge shape change of 10% at B o =0.002, g Q 
= 1.42x1 O' 5 g and C=0.001. The solid curve, the dotted curve and the 
dotted-dashed curve are the results for aspect ratios L/2R=2.6, 2.826 and 
3.024, respectively. From Zhang and Alexander [92], 

Fig. 14 Threshold intensity g* = g T R 4 /v 2 d T as a function of normalized frequency 

S = coR 2 /2v after Dewandre [57], Here g T is the maximum tolerable 
acceleration, d T the admissible displacement and r is the ratio of the drop 
density to that of the surrounding fluid. 

Fig. 15 Plot of effective diffusivity, Keff vs. Wormsely number W G , after Kurzweg 
[31]. 

Fig. 16 Dimensionless maximum velocity caused by impulse disturbances as a 
function of Fourier number F 0 = tic 2 /d, after Heiss et al. [64]. Here t is the 
time [s], k is the thermal diffusivity [cm 2 s’ 1 ], d is the diameter of the 
cylinder [cm]. The parameter Ra* is the Rayleigh number associated with the 
disturbance which is of dimensionless duration F 0 *= 5.5(1 0)* 3 . Ra* = (i T 
gL 3 AT/VK, where L is the cylinder length [cm], AT is the temperature 
difference between the ends of the cylinder [K], |3 T is the coefficient of 
thermal expansion [K* 1 ], g is the magnitude of the acceleration [cm s~ 2 ] and v 
is the kinematic viscosity [cm 2 s* 1 ]. 

Fig. 17 Maximum lateral solute segregation Ac/c« (where Coo is the far-field solute 

concentration in the melt) in the crystal as a function of the period of the 
gravitational oscillation, after McFadden and Coriell [66]. 



Fig. 18 


Fig. 19 


Fig.20 


Fig. 21 


Fig. 22 


Lateral non-uniformity in composition, £, plotted as a function of time for 
an oscillatory residual acceleration with a maximum magnitude of 3^2(10)' 6 
g and a frequency of 10' 3 Hz, acting parallel to the crystal-melt interface 
[67], The initial state was purely diffusive. 

Lateral non-uniformity in composition, £, plotted as a function of time for a 
residual acceleration consisting of a steady part with magnitude V2(10)" 6 g 
and an oscillatory part with a maximum magnitude of 3^2(10 )' 6 g and a 
frequency of 10' 3 Hz, acting parallel to the crystal- melt interface [67], The 
calculation was started from a steady flow associated with a V2(10)‘^ g 
acceleration acting parallel to the interface. 

Lateral non-uniformity in composition, plotted as a function of time for a 
multi-component disturbance consisting of a steady low g background plus 
two periodic components: g(t) = go + gi cos(27tlO~3t) + g 2 cos(27tl0"2t), 
where I g c I = V2(10)' 6 , I gi I = 3 V 2(10)' 6 and lg 2 I = 3^2(10 )' 5 [ 67 L 
The calculation was started from a steady flow associated with a V2(10)' 6 g 
acceleration acting perpendicular to the interface. 

Lateral non-uniformity 1^ I as a function of time for accelerations taken from 
a sample time series (see inset) constructed from data obtained on Spacelab 
3 [69]. The acceleration consists of a repeated “noise” segment (length 10 
seconds, 10' 1 Hz < /< 1 Hz) and a thruster firing event (length 10 seconds, 
10' 1 Hz </< 13 Hz). The latter is introduced at 10 and 80 seconds. 
Sensitivity of crystal growth from solution to various periodic disturbances. 
The solid line separates tolerable responses (10% or less growth rate 
fluctuations) from intolerable ones. For each case, the numbers correspond 
to the ratio g^AC/f v* , which is an estimate of the the ratio of bouyancy 
and inertial forces. After Nadarajah et al. [102]. 



Fig. 23 


Residual accelerations measured on orbit (in units of g = 9.8 m s' 2 ). (Data 
courtesy of Dr. H. Hamacher, DLR, and M.J.B. Rogers Univ. of Alabama 
in Huntsville.) The Spacelab 3 (SL-3) data are restricted to measured 
frequencies between 10'^ and 1 Hz. The points corresponding to Sled, 
Hop, Drop, Treadmill, Quiet, FPM ops, Stowage, refer to activities and 
experiments on the D-l mission. FPM ops stands for Fluid Physics Module 
operations, P RCS and V RCS refer to the primary and vernier thrusters. 
The drag and rotation entries correspond to accelerations arising from slow 
variations in atmospheric drag during an orbit and attitude changes 
involving rotation. 

Fig. 24 Acceleration sensitivity curves from selected examples discussed earlier 
superimposed on Fig. 23. Curves 1 and 3 correspond to the cases n=l, 
A=0.9999 and n=l, A =0.9 in Fig. 12 a, curve 2 corresponds to n=2 and 
A=0.9999 in Fig. 12 b. Curve 4 corresponds to the most sensitive case of 
zone shape change in Fig. 13, and curves 5 and 6 respectively correspond 
to the semi-conductor and metal melt growth experiments from Fig. 4. 
Curve 7 represents the thermo-diffusion experiment of Fig. 5. 



Amplitude Spectrum [10* 4 g] Residual Acceleration [10' 3 g] 





Amplitude Spectrum [1 0' 4 g] Residual Acceleration [g] 


0.012 

0.01 

0.008 

0.006 

0.004 

0.002 

0 



Time [sec] 



Frequency [Hz] 


Fig. 2 




Log Tolerable Acceleration [g] 




Log Tolerable Acceleration [g] 





Log Tolerable Acceleration [g] 





Grashof Number (Gr) 



Strouhal Number (St) 


Fig. 6 



Log Tolerable Acceleration [g] 



Log Frequency [Hz] 


Fig. 7a 




Log Frequency [Hz] 


Fig. 7 b 




log V g L a /D = log Pe g 


Fig. 8 





Fig. 9 


0.5 



Fig. 10 


L/2R 



Fig. 11 




Log Frequency [Hz] 


Fig. 12 a 



Log Tolerable Acceleration [g] 




ATION [g] 





Fig. 14 



Log g 





5000 



jxei/M 


Fig. 16 


Fourier Number 




Period 


Fig. 17 




8 

6 

4 



g. 18 






26 



g. 19 








Fig. 21 






[6] apniudiuv uopejaieoov 


Fig. 22 




N 





E 


CO 

Q_ 

O 

CD 

CD 

CO 

CO 

CO 


o 

V 

"O 

CD 

Q_ 

o 

Q_ 

O 

"O 

CO 

o> 

"cd 

Q_ 

£ 

o 

O 

cr 

o 

cr 

CD 

CO 

5 

o 

CO 

CO 


o 

h 1 

O 

U_ 

CO 

□_ 

> 

Q 

cr 

CO 

□ 

X 

• 

o 

◄ 

B 

+ 

< 

o 

o 

> 

• 



[6] uo!iej0|8oov Bon 


Fig. 23 


Log Frequency [Hz] 




r r • * 1 1 

[6] UO!iBJ0|9OOV 6on 


Fig. 24 


Log Frequency [Hz] 



APPENDIX 2. 


Accepted for publication. Physics of Fluids A: Fluid Dynamics, 1990 


Sensitivity of liquid bridges subject to axial residual 
acceleration 

Yiqiang Zhang and J. I wan D. Alexander 

Center for Microgravity and Materials Research 
University of Alabama In Huntsville 
Huntsville, Alabama 35899, USA 

It has become evident that the float zone crystal growth method 
and isothermal liquid bridges may be very sensitive to the residual 
acceleration environment of a spacelab. Indeed, the shape equilibria 
show a high degree of sensitivity and, thus, even the small steady 
acceleration associated with the effective low gravity environment of a 
spacecraft cannot be ignored. Using a slender-body approximation, 
the problem of determining the axisymmetric response of the shape of 
the free surface of a cylindrical liquid column bounded by two solid 
regions is modelled by 1-D system of non-linear equations. It is found 
that the sensitivity of the zone shape depends on the static Bond 
number, B 0 , aspect ratio and viscosity as well as the amplitude and 
frequency of the disturbance. The general trend is an increase in 
tolerable residual gravity with increasing frequency. At the 
eigenfrequencies of the zone, however, there are dramatic deviations 
from this trend. At these frequencies the tolerable residual gravity 
level can be two orders of magnitude lower at this frequency. For the 
cases considered the values of B 0 were taken to be 0.002 and 0.02 and 
the dimensionless viscosities, C = vfp/yRo) 1 / 2 , used were in the range 
0.001-0.01. Aspect ratios ranging from 96.3 to 82.7% of the Rayleigh 



limit were examined. For these cases, the frequencies associated with 
the lowest tolerable acceleration have been found in the lO^-lO' 1 Hz 
range. In terms of previously recorded and predicted residual 
accelerations the sensitive frequency ranges for the cases examined 
are 10' 2 -10’ 1 Hz and 1-10 Hz. Maximum tolerable residual gravity 
levels as low as 10' 6 g have been calculated. The effect of viscosity is 
seen to to increase the tolerable acceleration level for all frequencies. 
The equilibrium shape, as determined by the steady background 
acceleration, has a pronounced effect at low frequencies. A change in 
slenderness of the bridge markedly changes the sensitivity to residual 
acceleration as the Rayleigh limit is approached. 
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I. INTRODUCTION 

The low gravity environment of a spacecraft affords the 
opportunity to eliminate or minimize the often undesireable effects of 
gravity on fluid behavior. Both the fluid physics and the materials 
preparation community have an interest in the sensitivity of liquid 
zones to residual accelerations. A truly zero-gravity environment is 
never achieved and the residual acceleration which arises due to 
gravity gradient, atmospheric drag, crew motion, structural vibration 
etc. 1 ' 3 is characterized by an ever changing magnitude and orientation 
and affects the relative motion of any mass in the spacecraft reference 
frame. Whether a particular space flight experiment involves the study 
of isothermal zone shapes, or attempts to utilize the low gravity 
environment to maximize the length of the liquid zone in a float-zone 
crystal growth experiment, the effects of residual acceleration must be 
considered 4 - Experiments involving free surfaces of low viscosity 
liquids are sensitive to acceleration. Indeed, for some space 
experiments, unexpected results have been attributed to residual 
accelerations in several cases. 5 ' 9 Martinez 5 and Haynes 6 comment on 
the effects of an attitude change during the study of liquid column 
stability and liquid spreading kinetics on Spacelab-1. The residual 
acceleration environment also had a noticeable effect on stationary 
fluid masses. Roll rates of 0.13° per second were apparently "... 
sufficient to interfere strongly with some of the fluid configurations 
studied..." 8 Later experiments on D-l which involved static and 
rotating liquid columns also experienced some interference that has 


3 



been attributed to the residual acceleration environment. 7 In 
particular. Martinez recorded the following interference apparently 
caused by the residual acceleration environment. During the initial 
stages of the establishment of a long cylindrical column by liquid 
injection, a large amplitude deformation of the bridge was observed in 
response to shuttle maneuvers. Following this, a 10 cm liquid column 
exhibited a pronounced random oscillation. 

In an experiment seeking to excite a c-mode deformation 7 of 
the bridge an amphora mode 7 resulted. The residual acceleration 
environment seems to have caused a shift in the stability limits. The 
maximum residual acceleration recorded during these experiments 
was 10' 4 g. 

In Table I we have expressed the Bond number B n = (pg 0 Ro 2 
/y)gn/g 0 , where y is the surface tension, Ro is the radius of the ends 
of the zone, p is the liquid density and g n is the acceleration 
magnitude and g Q =980 cms* 2 , for the physical properties of liquid 
columns used in two of the space experiments discussed earlier. 
These estimates indicate that the experiment involving the long liquid 
bridges should be more sensitive. In contrast to Martinez, none of 
Padday's experiments were reported to have been affected by the 
residual accelerations. 

The numerical modelling of the response of a free liquid surface to 
time-dependent residual accelerations is complicated by the non- 
linearity introduced by the free boundary motion. Few studies have 
been conducted on this subject. Furthermore, the models to date are 
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based on the assumption of an inviscid bridge 1 ®* 11 , or have dealt with a 
linearized 1-D model. 12 Thus, we have formulated a simple non-linear 
1-D model of a viscous isothermal liquid zone and used it to examine 
the effect of g-jitter oriented parallel to the zone axis. This allows an 
examination of interaction between various components of a typically 
complex residual acceleration without the restriction that the 
response necessarily be linear. The model is based on a formulation 
developed by Lee 13 which has been successfully applied to several 
types of liquid bridge problems. 10.11.14-18 

We have examined the effects of single frequency g-Jitter 
disturbances. The zone response is described in terms of its deviation 
from its equilibrium shape. The equilibrium shape depends on the 
value of the steady background acceleration. The sensitivity is 
expressed in two ways: firstly zone breakage, defined in terms of the 
conditions under which the zone radius goes to zero, and secondly in 
terms of a shape change criterion which indicates when the zone 
deviates more than 5% from its equilibrium shape. An optimal 
searching scheme is used to find the sensitivity limits of the bridge. 

II. FORMULATION OF THE MODEL EQUATIONS 

Consider a slender liquid bridge (see Fig. 1) held between two 
parallel rigid circular disks of radius Ro separated by a distance L. The 

disks are aligned coaxially. The liquid is an isothermal incompressible 
Newtonian fluid with constant physical properties. It is held in place 
by capillary forces. Consistent with our aim to isolate the mechanical 


response of the surface to an acceleration g(t) acting along the 
cylinder axis we make the following additional assumptions: 

(i) Internal motion of the liquid bridge is caused only by 
capillary pressure gradients caused by deformation of the surface. 

(ii) The effect of the atmosphere around the bridge is 
negligible. 

(iii) The interaction of the residual acceleration environment 
with the experiment apparatus is such that the top and bottom disks 
vibrate in phase and the frame of reference for the calculations is then 
conveniently taken to be attached to the disks. 

The governing equations are rendered non-dimensional using Ro, 
(pRo 3 / y) 1 / 2 . and Ro(pR 0 3 /y)' 1/2 to respectively scale distance, time 
and velocity. Here p is the mass density, y is the surface tension. The 
dimensionless equations for mass and momentum transfer then take 
the following form: 

In the liquid (0 < r < R(z,t), -A < z < A), 

l3(ru) gw. 0 (!) 

r 3r 3z 


3u du 3u 3p , 

— + u— + w— = — + 

3t 3r dz dr 



1 3u d 2 u 
r dr + dz 2 



( 2 ) 


9w dw dw dp 

— — + u— — + w— — = - 
dt dv dz dz 


+ C 


id 2 W I dw 3 2 w\ 

+ raT + ^)- 


( 3 ) 


where p* is a reduced pressure given by p* = p + B 0 g(t)z, and C = 


6 


vfp/yR) 1 / 2 (with v the kinematic viscosity) is a measure of the relative 
strength of viscous and capillary forces. The function g(t) is given by 
g(t) = 1 + (Bi/B 0 )sincot. where co is the dimensionless circular 
frequency, B 0 = g 0 pR 0 2 /Y Is the usual static Bond number and Bi = 
gipRo 2 /yis a dynamic Bond number. We define the aspect ratio as A = 
L/2R 0 , where L is the zone length. The magnitudes of the steady and 
time-dependent acceleration components are given by g Q and gi 
respectively. The linear frequency (Hz) of the disturbance is given by / 
= 2rco)(pRo 3 /y) 1/2 - 

The boundary conditions at the surface r=R(z,t) are 


P - 


2Q_ 


f 3u , R 23w 
(l+R 2 ) or Z 3z 

[l±Bz_ - R zz ] + B 0 g(t)z. 


J3w 3u 
\9r 3zi 


(i + R ?r /2 r 


(4) 


02R z (|i 

\3r 


8w\ 
dz j 





0 . 


(5) 


which, respectively, represent the balance of normal and tangential 
force components at the surface. In addition, there is the kinematic 
condition 


3R 

at 


- u 


+ 



0. 


The conditions at z= ± A are 


( 6 ) 
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R(z,t) = 1 , w(r,z,t) = u(r,z,t) = 0, 


(7) 


which fixes the contact line between the disk and the liquid and 
ensures no slip of the fluid in contact with the disk. At r=0 the 
conditions 

u(O.z.t) = 0 , ^(O.z.t) = 0. (8) 

dr 

ensure axisymmetry. The formulation is complete upon specifying 
suitable initial conditions which we take in the form 

R(z,0)) =R<j(z), w(r.z,0)=u(r,z,0)=0, (9) 


along with the constraint that the volume remain constant; i.e. 



(z) dz 



(z,t) dz 


2 A. 


( 10 ) 


The above system of equations defines an unsteady free boundary 
problem since the location (and hence the shape) of the oscillating 
surface is a priori unknown and must be determined as part of the 
solution. 

In order to calculate the response of the shape of the liquid zone 
to time-dependent axial accelerations, we approximate the system of 
equations presented in the last section and formulate a one- 
dimensional model. This considerably reduces the complexity of the 
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problem and allows an examination of a relatively wide range of 
parameters without excessive computation times. Meseguer 10 used a 
one-dimensional model that is valid provided the slenderness A is 
large. (Recall that the slenderness of a static bridge is bounded from 
above by the Rayleigh limit i.e. A max = x.) If radial momentum effects 

are neglected 15 then eqs. (2) and (3) become decoupled and the 
following system of equations results 


l3(ru) 9w 
r dr + dz 


(ID 


dw 3w 9p* 

— + w — = - J — 

at dz dz 


+ c 


a*w 

dz T’ 


( 12 ) 


Since w is now depends only on z, eq. (11) is easily integrated and 
yields 


u(r,z,t) 


r 9w 
2 dz 


(13) 


Following Lee 13 and Meseguer 10 we take S=R 2 (z,t) and Q = 
S(z,t)w(z,t) and substitute (13) into the kinematic boundary condition 
(6). This yields 


as 

at 



(14) 
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which expresses conservation of mass at the surface. Similarly, eq. 
(12) can be recast in the form 


39 

a 

+ — 

(q 2 | . 

-s 9p * + sc f 

(Q\ 

(15) 

at 

dz 

‘ s 1 

dz dz 2 

\s) * 



which expresses the conservation of momentum in the axial direction 
for each cross section of the zone. At the surface, the component of 
force balance normal to the surface reduces to 


P = 


4S 




as \ 2 

dz I 


all* - “Hi 


+ 2S— 

dz 8z 2 *S 


4145 laz 


id - 


(2S 


jasf_ 

laz) 


2 ^2c 

S|§) + 

dz 2 


(16) 


where eq. (14) and the expressions for S and Q have been used to 
replace u and w. Together with the following boundary and initial 
conditions, eqs. (15) and (16) complete our approximate description 
of the physical system 

S(±A,t) =1, Q(±A,t) =0, 

S(±A,0) =S 0 (z), Q(±A,0) =0. (17) 

in. SOLUTION METHOD 

The time-dependent response of the shape of a liquid zone to 
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axial residual accelerations was calculated and curves describing the 
tolerable acceleration as a function of the frequency of the axial 
acceleration were obtained for a variety of conditions. The initial 
shape was determined from computed static shapes corresponding to 
the static Bond number Bo using a method similar to that employed by 
Meseguer and Sanz 15 . 

(a) Iteration procedure. In order to facilitate the solution of the 
system of eqs. (14-17) it is convenient to recast (15) and (16) in the 
form 


+!-(?)■ - s3 -f + cais) t2 (!) + cs (s, |- (I 1 • 

dt dz \ S / dz dz 2 \S> dz ‘S 


(18) 


P = 


4S 


+ 1—1 


2S as a 2 


(II + i 4S "if)' 1 ' 28 faff' S W ] + B ° 6(t)z - 


3S\ 2 \ dz dz 2 


where 


A (S) = S 


8S-(|S’ 2 

oz 

HiF 


( 20 ) 


and 


B (S) = S 


”§r-“fg 


( 21 ) 


(19) 


1 1 



The solution scheme is based on a fully implicit time stepping 
procedure. However, since the position of the surface is unknown an 
iterative procedure is required to obtain solutions to the differential 
equations and the location of the surface at each time. At the n+l th 
time step the iterations proceed as follows: 

(i) Guess the surface shape for the initial iterate (m=l). 

(ii) The pressure p* at the m+1 th iterate is calculated in the 
following way. All terms in (19) are calculated from the previous (m th ) 
iterate except for the leading viscous term. This is calculated from 
the solutions obtained at the previous (n th , n-l th ) time steps using an 
Adams-Bashforth two-step approximation. 

(iii) Q is then obtained upon solution of the momentum equation 
(18) and boundary conditions (17) using a Briley- M c Donald implicit 
method 19 ' 21 . 

(iv) The surface shape at the m+l th iterate is then obtained from 
a time-centered discretization of (14). 

(v) The surface shape is then compared with the initial guess (or 
the shape calculated at the m th iterate). The procedure is repeated 
until the condition 


maxl- 


isr 1 

Sf 


sr 1 


m 


< e 


( 22 ) 


where e «1, is satisfied. Here a subscript 'i' denotes the value of s at 
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the i th point, while the superscript denotes the iterative step. The 
above procedure is repeated at each time step. 

The Briley-M c Donald method used to solve (18) involves the 
evaluation of the non-linear advection term at the n+l th time step 
using 


aioVs) 

n+1 

afoVs)' 

n 

+ At 

i 

a (a(9 2 /s)\' 

dz ■ 

i 

dz 

dz\ at I 


in 


O ((At) 2 ) 


(23) 


where the time derivative has been approximated using a forward 
difference and the space derivative with a centered difference. A 
Crank-Nicolson scheme is used to increase the accuracy of the 
solution. The discretized form of (18) then becomes 




+ B (S)|- 
dz 




A(s) |r(- 


+ B( S)|- 

az 



(24) 


where a centered difference is used for spatial discretization. Eq. (24) 
is unconditionally stable and has second order accuracy with a 
truncation error of 0[(Ax) 2 , (At) 2 ]. 

In addition to the boundary conditions (17) we also need to 
evaluate the pressure at the disks. This is achieved as follows. First we 
use the fact that ds/dt = 0 at z=± A which implies dQ/dz = 0 (from (14)) 
and that dQ/dt =0 at z =± A. We can then obtain the pressure from (18) 
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evaluated at z =± A. This requires the evaluation of dp/dz and 3 2 Q/9z 2 
a t z =- A. The pressure derivative is approximated using the following 
standard three-point forward difference expression: 


(&) 

\3z /2 = +^\ 


= ± 


P2 - 4px + 3pp 
2Az 


+ O (Az3) 


(25) 


where, po - p(±A), pi = p(± A + Az), p2 = p(± A + 2Az) and Az is the 
spatial step size. Similarly, at z= ± A, d^Q/dz 2 is evaluated using 


« 5z2 lz=±A 


8Q1 - Q 2 - 7Q 0 
2Az 2 


+ O (Az 4 ). 


(26) 


where Q 1 - Q(± A + Az), Q 2 = Q(± A + 2Az) and the fact that dQ/dz = 0 
has been employed. Upon substitution of (25) and (26) into (18) and 
rearranging terms we obtain p 0 with second order accuracy. 

The system of discretized equations and boundary conditions is 
solved using a Thomas algorithm (method of factorization) 19 . 

(b) Initial iterate For this method of solution the choice of 
surface shape for the initial iterate at each time step affects the speed 
of convergence. For a poor choice the iteration may even fail to 
converge. To ensure that the initial iterate is reasonably close to the 
solution, the following formula is employed 



14 


For the computations presented here the number of iterations 
required was six or less whenever (22) was satisfied for the initial 

iterate. 

( c ) Optimal searching scheme In order to define the sensitivity 
limits for the liquid zone shape an optimal searching scheme was used 
to delineate the boundary between the regions in parameter space for 
which solutions either do or do not satisfy our sensitivity criterion. 
The scheme is described in the flow chart shown in Fig. 2. Eight to 
ten calculations were typically needed to obtain each point on the 

sensitivity curves. 

(d) Accuracy of the numerical scheme An investigation into the 
accuracy of the results revealed that the number of timesteps per 
period of the disturbance should be at least 160. Increasing the 
number of steps to 400 per period did not, however, result in any 
significant change in our results even at high frequencies (see Table 

II). 

The spatial accuracy also proved to be an important 
consideration. Results obtained using the MacDonald-Briley scheme 
required 97 spatial points. 

In order to test the adequacy of the 1-D model we compared results 
from 1-D calculations to results obtained by Sanz 16 for breaking times 
of an inviscid liquid bridge subject to an initail deformation of the 

surface area S given by 
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S(z.O) = 1 + e(e-2)sin|S2-J, Q(z,0) = 0. 

The results of the comparlsopn are given in Table III. Our calculations 
are in reasonable agreement with those of Sanz. Table IV lists a 
comparison between our results for the first eigenfrequency, those 
predicted by the linearized analysis of Meseguer 11 and experimentally 
obtained values 11 . Our calculated values lie between the Inearized 
results and the experimentally obtained values. Finally, we note that 
the analysis of Langbein 12 , employs an expression for a damped linear 
oscillator together with natural frequencies of infinite liquid columns 
obtained by Bauer 22 and will overestimate the natural frequencies of 
finite bridges. This has been well documented by Sanz 16 in a thorough 

comparison of various methods used to simulate liquid bridge 
phenomenon. 

IV. RESULTS AND DISCUSSION 

The calculations were undertaken for a disk radius of 0.0175 m 
and liquids with the following physical properties; p = 920 kgnr 6 . y 
=0.02 Nm* 1 . and kinematic viscosities in the range 6xl0* 7 - 6x1 O' 5 
m 2 s* 1 . These viscosities cover a range of typical experimental 
materials (water, silicone oils etc., used in experiments 10 ). The 
sensitivity criteria chosen to characterize the response of the bridge 
represent two extremes. The first is the deviation of the the bridge 
shape by more than 5% from its equilibrium radius, the second is 
breakage of the bridge. The effects of viscosity, background steady 
acceleration (static Bond number, B 0 ) and slenderness (A) were 
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examined and the results are presented below. Each point on these 
curves was obtained after running the calculation for times 
corresponding to more than 10 periods of the driving force. 

Figures 3 and 4 show the effect of viscosity on the tolerable 
acceleration. g t . for a frequency of 5 Hz at fixed A and B 0 . for both 
sensitivity criteria. The tolerance increases as the viscosity is 
increased. Notice, however, that for values of the dimensionless 
viscosity parameter C less than 10‘* the increase in tolerance is slight 
even though the viscosity has been increased by two orders of 
magnitude. Only when C approaches 10* 1 does the increase in 
tolerance become significant. Figures 5 and 6 further illustrate this 
point. 5 Hz residual accelerations greater than 0.03 g would not 
normally be expected to occur in spacecraft laboratories so we have 

not examined cases with C> 0. 1 . 

The value of the static Bond number has a more pronounced 

effect on sensitivity as seen in Figs. 7 and 8. This is limited, however 
to the lower frequency range. With increasing frequency the 
differences are less pronounced. This is particularly so for the shape 
change criterion. It is interesting to note the additional deviations 
from the general trend for the case of zone breakage between 10' 1 and 

1 Hz. 

A small decrease in the slenderness A changes the magnitude of 
tolerable acceleration markedly when A is close to A ma x- This is 
illustrated in Figs. 9 and 10. The change in tolerable acceleration 
magnitude becomes less significant for values of A less than 80% of 
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the critical value. 

Of most practical interest is the sensitivity of a given zone (i.e. A. 
B 0 , and C fixed) as a function of the frequency of the disturbance. 
This is shown in Figs. 5-10. For both criteria the general trend is an 
increase in tolerable residual gravity with increasing frequency. At the 
eigenfrequencies of the zone there are, however, dramatic deviations 
from this trend. Here the liquid zone is extremely sensitive in 
comparison to both higher and lower frequencies. The tolerable 
residual gravity level can be two orders of magnitude lower at this 
frequency. For the values of A and B 0 (g 0 ) considered these 
frequencies are found in the neighborhood of 10' 1 Hz. Associated 
maximum tolerable residual gravity levels as low as 10' 3 to 10' 6 g have 
been calculated. In addition, at frequencies higher than the most 
sensitive one there are less dramatic deviations from the general 
trend. The frequencies at which these deviations occur are the 
eigenfrequencies for the zone. As the frequency of the forcing 
function approaches an eigenfrequency the zone becomes more 
sensitive. The lowest tolerable acceleration occurs in conjunction 
with the smallest eigenfrequency. We note that owing to our 
assumption that the top and bottom disks vibrate in phase, the 
eigenfrequencies locally exhibiting the most sensitive response are 
those associated with eigenmodes having an odd number of nodes (i.e. 
an even number of half waves). This is in agreement with the analysis 
of Mesegeuer 11 . It is interesting to compare our results with linear 
oscillator model of Langbeini^ , Table V gives a comparison between 
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the linear oscillator model and our 1-D model. While the results are 
qualitatively similar there are three distinct differences. Firstly, for the 
case considered the linear oscillator model overestimates the 
sensitivities by almost two orders of magnitude. Secondly, as 
mentioned previously, the eigenfrequencies predicted by our model 
are smaller. The third difference is that Langbein’s model implicitly 
assumes that the axial vibration will excite deformations with an both 
odd and even number of nodes (this is precluded in our model). If the 
top and bottom disks do vibrate out of phase, then our results will 
underestimate the response of the liquid bridge at eigenfrequencies 
corresponding to eigenmodes with an even number of nodes. 

Table VI lists the eigenfrequencies (which can be determined 
independently 11 ) and the associated maximum tolerable acceleration, 
gt, for selected inviscid cases. For the invlscid zones, we found that 
whenever the higher eigenfrequency is close to an integer multiple of 
a lower eigenfrequency the sensitivity of the zone to disturbances with 
the higher frequency is increased. This noticeable increase in 
sensitivity is caused by non-linear interaction and excitation of the 
lower mode. 

The spatial deformation of the liquid zone also varies as a 
function of the forcing frequency. Figures 11 and 12 illustrate the 
evolution of the zone shape as a function of time for selected cases. 
Note that, owing to the small value of B 0 . for these cases the initial 
surface shape is visually indistinguishable from the surface of a right 
circular cylinder. While higher forcing frequencies are associated with 
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more nodes, the shape of the zone is complicated by the excitations of 
the fundamental mode and by non-linear interactions. 

In summary, our results indicate that the zone is most sensitive 
to accelerations with frequencies close to or equal to the lowest 
natural frequency of the zone. It is useful to assess the sensitivity of 
the zone in terms of predicted space-station and/or spacelab 
environments and with previous acceleration measurements. For the 
liquid zone cases examined we have seen that frequencies in the 10' 2 - 
10’ 1 Hz range appear to be the most sensitive. Figure 13 shows 
sensitivity curves selected from Figs. 9 and 10 together with 
acceleration as a function of frequency taken data recorded on orbit. 
The amplitudes of low frequency (< 10' 2 Hz) accelerations predicted 
for the space station should not exceed levels of 10‘ 6 g 23 . Higher 
frequencies can be associated with acceleration magnitudes of up to 
10-2 g24. i n terms of these predicted levels, or those measured on 
past missions 3 - 25 , the practical sensitivity range is restricted to 
disturbances with frequencies in the range 10* 2 - lCH Hz and 1-10 Hz. 
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TABLE I. Dynamic Bond numbers estimated for two D-l liquid zone 
experiments in which 5 cs silicone oil was used. 



Martinez 5 

Padday 6 

p [g cm' 3 ] 

0.92 

0.92 

L [cm] 

10 

1.5 

c [dyne cm -1 ] 

0.2 

0.19 

Bn ' 

4.5(10)' 4 gn/gO 

7.1(1 0) -3 gn/gO 
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TABLE II. Effect of changing At on the tolerable acceleration 


gt- 


At 

CO 

A 

Bo 

gt 

0.00125 

0.0005 

31.4 

31.4 

2.6 

2.6 

0.002 

0.002 

0.255 

0.256 
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Table III. Comparison of breaking times (T b ) and partial breaking 
volume (V b ) with those of Sanzie ( Tb *. y b *). A is the slenderness, e is 
the initial perturbation amplitude and Az is the spatial step size used 
in our calculations. 


A 

e 

Az 



v b 

v b * 

3.07 

0.2 

0.128 

21.02 

21.98 

0.8525 

0.8554 

3.14 

0.2 

0.130 

13.71 

13.62 

0.8474 

0.8574 

3.20 

0.2 

0.134 

11.82 

11.52 

0.8427 

0.8442 

3.10 

0.3 

0.129 

7.97 

8.4 

0.8525 

0.8539 

3.10 

0.4 

0.129 

4.41 

4.84 

0.85945 

0.8566 
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Table IV. Comparison of the dimensionless first eigenfrequency «o 
computed using our method, the linearized method of 
Meseguer 1 and experimental results 1 1 (coo e ) 


A 

COo 

(0*o 

Woe 

2.571 

0.312 

0.348 

0.30 

2.714 

0.248 

0.259 

0.23 

2.854 

0.174 

0.222 

0.19 
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Table V. Most sensitive frequencies and associated tolerable 
acceleration [Hz] calculated using the 1-D [f 0 , g tol ) model and 
Langbein’s [f q g*tol) linear oscillator model 12 . For the 1-D model Bo = 
0.002 in all cases. Langbein's model is independent of Bo. 


A 

/O [Hz] 

f o[Hz] 

gtol [gl 

g*tol [g] 

v [m 2 s' 1 ] 

3.024 

0.0334 

0.0663 

6.32xl0* 7 

1.11x10-9 

5xl0- 7 

2.826 

0.0070 

0.1224 

1.16x10-6 

2.18x10-9 

5xl0- 7 

2.6 

0.1049 

0.1857 

3.26x10-6 

3.6x10-9 

5x10-7 

2.6 

0.1049 

0.1857 

3.53x10-6 

3.6xl0- 8 

5x10-6 
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TABLE VI. Tolerable acceleration (gt) and associated eigenfrequencies 
(co*) for two of the cases (shape change criterion). 


A = 3,024.Bn = 0.002. C =0 

m* 0.21 4.0 

11.0 

20.5 

33.5 

gt (m s' 2 ) 6.0xl0* 6 6.09xl0' 4 

1.88xl0- 3 

5.10x10-3 

5.06x10-3 

A = 2.826.B^= (LQ02. C =0 

0)* 0.44 4.7 

12.56 

24 

38.5 

gt (m s’ 2 ) 1.09x10-5 7.79x10*4 

1.96x10-3 

4.76x10-3 

5.78x10-3 
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FIGURE CAPTIONS 


FIG. 1 Idealized liquid zone configuration. 

FIG. 2 Computational flow chart. 

FIG. 3 Dimensionless viscosity parameter C vs. tolerable 

acceleration for breakage of the bridge at f =5 Hz, A=2.6, Bo=0.002 
and go=1.42xl0’5 g. 

FIG. 4 Effect of viscosity on tolerable acceleration for the shape 
change criterion at f=5 Hz. A=2.6, Bo=0.002 and go=l. 42x10*5 g. 

FIG. 5 Curves of tolerable acceleration versus frequency for 
breakage of the bridge at Bo=0.002, g Q =l. 42x10-5 g and A=2.6. 

The solid curve and the dashed curve are the results for C =0.001 
and 0.01, respectively. 

FIG. 6 Curves of tolerable acceleration versus frequency for shape 
change criterion at Bo=0.002, g 0 = 1.42x10-5 g and A=2.6.The solid 

curve and the dashed curve are the results for C=0.001 and 0.01, 
respectively 

FIG. 7 Curves of tolerable acceleration versus frequency for shape 
change criterion at A=2.6 and C=0.001. The solid curve and the 
dashed curve are the results for Bo=0.002 (g Q = 1.42x10*5 g) and 
Bo=0.02 (go= 1.42x10*4 g), respectively. 

FIG. 8 Curves of tolerable acceleration versus frequency for 
breakage of the bridge at A=2.6 and C=0.001. The solid curve and 
the dashed curve are the results for Bo=0.002 (g 0 = 1.42x1 O' 5 g) 
and Bo=0.02 (g 0 =l. 42x10*4 g), respectively. 


FIG. 9 Curves of tolerable acceleration versus frequency for 
breakage of the bridge at Bo=0.002, go= 1.42x1 0*5 g and C=0.001 

The solid curve, the dotted curve and the dotted-dashed curve are 
the results for A=2.6, 2.826 and 3.024, respectively. 

FIG. 10 Curves of tolerable acceleration versus frequency for the 
shape change criterion at Bo=0.002, g 0 =1.42xlO* 5 g and C=0.001 

The solid curve, the dotted curve and the dotted-dashed curve are 
the results for A=2.6, 2.826 and 3.024, respectively. 

FICj. 1 1 Shape change as a function of time for a liquid zone with 
A=2.6, C=0.001 and g(t) =1 + (Bi/B 0 )sin(cot). The dimensional 
frequency, J , is 0.06Hz, and B^/Bq =12. For the physical properties 

considered, the maximum oscillating acceleration magnitude is 
1.8xl0- 4 g. 

FIG. 12 Shape change as a function of time for a liquid zone with 
C=0.001 and g(t) =1 + (Bi/B 0 )sin(cot). The dimensional 
frequency, f , is 3 Hz, and Bi/B 0 =916. For the physical properties 

considered, the maximum oscillating acceleration magnitude is 
1.3x10-2 g. 


FIG. 13 Residual accelerations measured on orbit (in units of g = 
9.8 m s’ 2 ) and selected sensitivity curves. Curves 1 and 2 
correspond to A= 2.826 and 3.024 from Fig. 10. Curve 3 
correspond s to A= 3.024 from Fig. 9. Acceleration data courtesy of 
Hamacher, DLR, and M.J.B. Rogers Univ. of Alabama in 
Huntsville.) The Spacelab 3 (SL-3) data are restricted to measured 
frequencies between 10’ 3 and 1 Hz. The points corresponding to 
Sled. Hop, Drop. Treadmill, Quiet, FPM ops (Fluid Physics Module 
operations). Stowage, refer to activities and experiments on the D- 
1 mission 23 . P RCS and V RCS refer to the primary and vernier 
thrusters. The drag and rotation entries correspond to 
accelerations arising from slow variations in atmospheric drag 
during an orbit and attitude changes involving rotation. 
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